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Abstract 

Our aim is to give a self-contained review of recent advances in the analytic de- 
scription of the deconfinement transition and determination of the deconfinement 
temperature in lattice QCD at large A^. We also include some new results, as for 
instance in the comparison of the analytic results with Montecarlo simulations. We 
first review the general set-up of finite temperature lattice gauge theories, using 
asymmetric lattices, and develop a consistent perturbative expansion in the cou- 
pling (3s of the space-like plaquettes. We study in detail the effective models for 
the Polyakov loop obtained, in the zeroth order approximation in /3s, both from 
the Wilson action (symmetric lattice) and from the heat kernel action (completely 
asymmetric lattice). The distinctive feature of the heat kernel model is its relation 
with two-dimensional QCD on a cylinder; the Wilson model, on the other hand, 
can be exactly reduced to a twisted one-plaquette model via a procedure of the 
Eguchi-Kawai type. In the weak coupling regime both models can be related to 
exactly solvable Kazakov-Migdal matrix models. The instability of the weak cou- 
pling solution is due in both cases to a condensation of instantons; in the heat kernel 
case, this is directly related to the Douglas-Kazakov transition of QCD2. A detailed 
analysis of these results provides rather accurate predictions of the deconfinement 
temperature. In spite of the zeroth order approximation they are in good agreement 
with the Montecarlo simulations in 2 + 1 dimensions, while in 3 + 1 dimensions they 
only agree with the Montecarlo results away from the continuum limit. 
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1 Introduction 



In the last twenty years the lattice regularization has proved to be a very power- 
ful tool to understand and describe the non-perturbative features of non-abelian 
gauge theories. However, while impressive results have been obtained by means 
of Montecarlo simulations, very few progresses have been achieved with analytic 
techniques. Lattice Gauge Theories (LGT in the following) can be solved exactly 
in two dimensions for any gauge group, but become unaffordably complex in more 
than two dimensions, even in absence of quarks. Moreover, most of the approxi- 
mation techniques which are usually successful in dealing with simpler statistical 
mechanical systems, like (suitably improved) mean field methods or strong coupling 
expansions turn out to be less useful in the case of LGT. A remarkable exception 
to this state of art is represented by the large N approximation which is able to 
keep the whole complexity of the finite N models. Unfortunately, even in the large 

limit (despite the fact that, as we shall discuss below, some major simplifications 
occur) it is not possible to give exact solutions (the so called "Master Field" ) to the 
Lattice SU(A^) models. Notwithstanding this, several interesting results have been 
obtained in the past years even without the explicit knowledge of the Master Field. 

In this review we shall deal with one of the most interesting features of non 
abelian gauge theories: the presence of a deconfinement transition at finite tem- 
perature. We shall apply large techniques to finite temperature LGT, with the 
aim of constructing the phase diagram of the model, locate the critical points and 
identify their order. We shall deal for most part of the review with the pure gauge 
theory (namely without quarks); only in section ^.2| we shall comment on the phase 
diagram in presence of quarks. We shall keep d, the number of space-like dimensions 
as a free parameter, and shall study in particular the cases d = 2 and d = 3 which 
are the most interesting from a phenomenological point of view and for which there 
exist Montecarlo simulations to compare with our predictions. 

This paper is an update rather than a complete review on this topic. We shall 
mainly concentrate on the most recent results (say, of the last five years). There 
are some very good reviews both on lattice gauge theories and on large A^ models 
which cover most of the results obtained up to the second part of the eighties. In 
particular, for the large A^ approach we suggest the reviews of S. Coleman 0, S.R. 
Das and the recent contribution of Campostrini, Rossi and Vicari For a 
throughout introduction to Lattice Gauge Theories, with a good discussion both of 
the large A^ limit and of the finite temperature regularization, we refer to 

During the eighties most of the efforts in the study of finite temperature large A^ 
LGT where devoted in the following two directions: a) Eguchi-Kawai (EK) type ^ 
models; b) dimensionally reduced model. Both these approaches have advantages 
and drawbacks. Twisted or Quenched EK models maintain the whole complexity 
of the theory and reduce it to models of just {d+ 1) independent matrices. However 
these models could not be solved exactly, and the only fruitful approach has been so 
far to extract numerical results by using Montecarlo simulations. We shall discuss 
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below some of these results. Let us mention here that the reliability of Montecarlo 
simulations based on the TEK model is strongly affected by the presence of a 
first order bulk phase transition which shadows the true deconfinement transition. 
Finding a precise characterization of the TEK or QEK phase diagram, from the 
results of the Montecarlo simulations only, is a difficult and open problem. 

On the other side, by doing some rather crude approximation, it has been pos- 
sible to construct dimensionally reduced matrix models which are exactly solvable. 
The price to pay in this case is that the simplifications needed to reach the exact 
solvability are so strong that the resulting models show sometimes only a very faint 
similarity with the original models. 

In the last few years it has been realized that some insight into the large 
structure of LGT can be obtained by using the results that have been accumulating 
in the meantime on the large solution of a variety of matrix models. These models 
have found important applications in various contexts. Let us mention among the 
others: the random matrix description of various physical systems (ranging from 
quantum wires to chaotic systems), the exact solution of two dimensional LGT's 
P, ^ |T0|, [Tl|] and the related string like behaviour , the large N matrix approach 
to two dimensional quantum gravity [13| and the "induced QCD" models proposal 
by Kazakov and Migdal il^ 



The present review is mostly concerned with the results ^7\, |T8[ obtained 



in the last few years in finite temperature LGT by using the matrix models tech- 
niques mentioned above. We shall show that by using both some recent results 
of two dimensional QCD and the exact solution |]TU| of the induced QCD mod- 
els in any dimension, one can solve dimensionally reduced models (of the type b) 
listed above) that involve far less crude approximations than those discussed in the 
eighties. By using these exact results we can also solve, within a very good approx- 
imation, a simplified TEK model and understand its phase diagram. Our results 
can be compared with those obtained with Montecarlo simulations in the N = 2 
and N = 3 cases. There are by now many rather precise numerical estimates of the 
deconfinement temperature for these models both in d = 2 and d = 3. We shall 
show that our results are rather good if compared with the simulations at ci = 2, 
while they are in general poor (except for the lowest values of rit, the lattice size in 
the time direction) in d = 3. We shall give some arguments to explain this failure 
and indicate how our results could be improved. Let us stress however that our 
aim is not to compete with the Montecarlo simulations, which certainly remain the 
most powerful tool to obtain quantitative result in LGT. Our idea is rather that any 
new analytic result in this context can teach us quite a lot on the non-perturbative 
regime of non abelian gauge theories. Moreover, even if our results are often only 
crude approximations of the exact ones, the pattern of the phase diagram that we 
obtain has good chances to be the correct one. Finally let us stress that there are 
situations in which the methods discussed here could even become competitive with 
the Montecarlo simulations. This is the case, for instance, when the model contains 
more that one coupling constant (like the mixed fundamental-adjoint action), and 
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the phase diagram is so complex that it is difficult to study it numerically. 

The paper is organised as follows. For the sake of being as self-contained as 
possible, we devote a rather large introductory section (sec. 0) to describe the gen- 
eral setting of finite temperature LGT in a form which is suitable for our purposes. 
In sec. ^ we review in some detail the derivation of refined effective action for the 
Polyakov loops, both in the case of a LGT described by the heat kernel action and 
by the Wilson action. In sec. | the heat kernel effective model is discussed; its 
phase diagram is described, and a numerical estimate for the deconfinement tem- 
perature is obtained. In sec. |^ the phase diagram is re-discussed as a function of 
the space dimensionality and in the case in which the gauge fields are coupled to 
certain external static sources or to an external "magnetic" field. Sec. |] is devoted 
to the analysis of the effective model obtained from the Wilson action; in partic- 
ular we discuss its twisted reduction a la Eguchi-Kawai, that gives rise to a very 
interesting one-plaquette matrix model. In sec. |^ the results for the deconfinement 
temperature are discussed and compared with the available Montecarlo simulations. 



2 General Setting 
2.1 Asymmetric Lattices 

Let us consider a finite temperature lattice gauge theory (LGT) with gauge group 
SU(A^), defined on a c? -|- 1 dimensional cubic lattice. In order to describe a finite 
temperature theory we require that one dimension (which we shall call "time" from 
now on) is compactified with periodic boundary conditions and compactification 
length j;, if T is the temperature. We shall assume that in the other "space" 
dimensions the extension of the lattice is infinite, or anyway much greater than j;. 
If we denote by Nt the number of lattice spacings in the time direction, and at the 
corresponding lattice spacing, then we have j; = Ntat- It is convenient to work with 
an asymmetric lattice, namely with different lattice spacing in the time and space 
directions. Let us denote by the lattice spacing in the space directions. The ratio 
between at and defines the asymmetry parameter p: 

p=- = TNta, . (1) 
at 

We shall require that different values of p correspond to different, but equivalent, 
lattice regularization of the same model. In order to implement such requirement 
we introduce different bare couplings in the time and space directions. Let us call 
them j3t and (3s respectively. The Wilson action is then 

= iV' E f E AG'o.(x, t) + E ^G^,{x, t) J , (2) 

x,t \ i i<j ) 



6 



where here and in what follows we denote by Got and Gij the time-like and space-like 
plaquette variables: 

Goi{x,t) = V{x,t)Ui{x,t+l)V\x + %t)Ul{x,t) , 

Gij{x,t) = Ui{x,t)Uj{x + i,t)Uj{x+],t)U}{x,t) (3) 

and by GQi{x,t) and Gij{x,t) the real part of the suitably normalized traces: 

Goi{x,t) = ^ReTrG'oi(x,t) , Gij{x,t) = ^ReTTGij{x,t) . (4) 

In eq. we have denoted by Ui{x,t) the link variables in the space directions 
{i G {1,2, d}) and by V{x, t) the link variables in the time direction. The com- 
ponents of X and t are integers, with t periodic modulo Nt. The normalization of 
the couplings Ps and Pt has been chosen, by extracting a A^^ factor in front of S'w, 
in such a way to have a smooth large N limit. 

For a given p the relation between the bare couplings Pt and Ps can be obtained 
by requiring that the Wilson action reproduces in the naive continuum limit 
a gauge theory with the same coupling constant g for all components of the field 
strength. This leads to the following equations, which relate p and the bare gauge 
couphng g to Pt and Ps- 



a'-yPsPt, P=\l^- (5) 



Ng' 



It is clear from eqs. ([1|) and (^) that equivalent regularizations with different 
values of p require different values of Nf. Hence, to maintain the equivalence, Nt 
must be a function of p according to eq. 

Among all these equivalent regularizations a particular role is played by the 
symmetric one, which is defined by: 

(from now on we shall distinguish the symmetric regularization from the asymmetric 
ones by eliminating the subscripts t and s in P and a). By comparing eq.s 
we see that all regularizations are equivalent to a symmetric one provided 

P = pPs = -, (7) 
p 

Nt{p)=pnt, (8) 

where nt is the number of links in the time direction with a symmetric regularization: 
nt = Nt{p=l). 

Notice however that these equivalence relations have been derived in the naive or 
"classical" continuum limit, and quantum corrections are in general present. These 
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corrections were studied and calculated in the (3+1) dimensional case by F.Karsch 
in pO|. They lead to the following expressions: 



l3t = p{l3 + cAp)) , (9) 

A = ^i^^ . (10) 

P 

The quantum effects are encoded in the two functions Co-(p) and c,-(p)[|. We shall 
be interested in their behaviour at large p, that can be extracted from |^ and is 
given by the asymptotic expansion: 

c.,r^a'+^ + ... . (11) 



The numerical values of the a's in the large limit can be obtained from |20 
and they are given by: a° = -0.2629; a]. = 1/4; a° = 1/4; = 1/4. 



2.2 Center symmetry and the Polyakov loop 

The major consequence of the periodic boundary conditions in the time direction is 
the appearance of a new global symmetry of the action, with symmetry group the 
center C of the gauge group (in our case Z^). 

This symmetry can be realized as follows. Let us transform all the timelike links 
of a given space-like slice with the same element Wq belonging to the center of the 
gauge group. 

V{x,t) ^WoV{x,t) Vf, t fixed. (12) 

The space-like plaquettes are not affected by the transformation, while in each 
timelike plaquette two contribution appear: Wq and Wq^. Since they belong to the 
center, they commute with all other matrices in the plaquette and can be moved so 
as to cancel each other. So the Wilson action is invariant under such transformation. 
The important point is that, due to the periodic boundary conditions in the time 
direction, it is impossible to re-absorb this global twist by means of local gauge 
transformations. . 

A second consequence of the periodic boundary conditions is that it is possible to 
define gauge invariant observables which are topologically non-trivial. The simplest 
choice is the Polyakov loop, defined in terms of link variables as: 

Nt 

P{x) =Til[V{x,t) . (13) 
t=i 

^Lct us notice, to avoid confusion, that our functions Cg- and Cr correspond to those of [ pO{ 
multiphed by the factor 2/N, which ensures a smooth hmit as N oo. 
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In the following we shall often use the untraced quantity P{x), defined as: 



P{x)^Y[V{x,t) , 



(14) 



which will be referred to as "Polyakov line" . 

The relevant feature of the Polyakov loop is that it transforms under the above 
discussed symmetry as follows: 



thus it is a natural order parameter for this symmetry. It will acquire a non zero 
expectation value if the center symmetry is spontaneously broken. In a pure LGT 
the Polyakov loop has a deep physical interpretation, since its expectation value is 
related to the free energy of a single isolated quark. 

Hence the fact that the Polyakov loop acquires a non-zero expectation value 
can be considered as a signature of deconfinement and the phase transition which 
separates the regime in which the center symmetry is unbroken from the broken 
symmetric phase will be the deconfinement transition. 

General arguments show that, in > 1, finite temperature gauge theories admit 
such a deconfinement transition for some critical value of the temperature T = 
Tc, separating the high temperature, deconfined, phase (T > Tc) from the low 
temperature, confining domain (T < Tc). In the following we shall be interested in 
the phase diagram of the model as a function of T, and we shall make some attempt 
to locate the critical point T^. 

A peculiar feature of the behaviour of the Polyakov loop in the high temperature 
phase is that its (non- zero) expectation value is not a generic element of SU(A^) but 
tends to fluctuate around one of the elements of the center. These fluctuations 
become smaller and smaller as the temperature increases and finally in the infinite 
temperature limit the expectation value of the Polyakov loop exactly becomes an 
element of C (see for instance |^ for a discussion of this point). 

This feature will be relevant in the following, when we shall describe these small 
fluctuations of the Polyakov loop in the high temperature phase by using a suitable 
generalization of the Kazakov-Migdal model. 

2.3 Svetitsky— YafFe conjecture 

The peculiar role played by the Polyakov loops in the above discussion, suggests to 
use some kind of effective action for the Polyakov loops to study the deconfinement 
transition and, more generally, the physics of finite temperature LGT. 

Let us make this statement more precise. Constructing an exact effective action 
for the Polyakov loops, equivalent to the original LGT is clearly impossible (even in 
the simplest possible, non-trivial, LGT, namely the (2+1) dimensional Ising gauge 
model) since it would require to integrate out exactly all the space-hke gauge degrees 



P(f) WoP{x) ; 



(15) 
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of freedom of the original model. However one can try to approximate somehow 
this integration and, at the same time, treat exactly the timelike degrees of freedom 
which are related to the Polyakov loops. The discussion of sec. 2.2 then tells us that 
this type of approximation is the one which better preserves the finite temperature 
behaviour of the original model. The resulting approximate effective action will keep 
all the symmetries of the original model relevant to the deconfinement dynamics 
and as a consequence it will give a faithful qualitative description of this transition. 
We can also hope to have a quantitative agreement between the results obtained 
with the effective action and those of the gauge model. This agreement will become 
better and better as we improve our approximation in the space-like degrees of 
freedom. 



This approach was proposed and discussed in |21|] and accordingly we shall re- 
fer to it, in the following, as the "Svetitsky-Yaffe program". The most appealing 
feature of this program is that, by resorting only to some general results of sta- 
tistical mechanics, before doing any calculation, one can obtain several interesting 
predictions on the phase diagram of the model. Let us see some of them: 

a) If the original gauge theory lives in (d+l) dimension, then the effective theory 
for the Polyakov loops is a ci-dimensional spin system with symmetry group 
the center C of the original gauge group. 

b ) The deconfinement transition of the original gauge model becomes the order- 
disorder transition of the effective spin system. The ordered phase of the spin 
model corresponds to the deconfined phase in the original gauge theory. In 
this phase the Polyakov loop acquires a non-zero expectation value. 

a) This effective theory would obviously have very complicate interactions, but 
Svetitsky and Yaffe were able to argue that all these interactions should be 
short ranged. As a consequence, if the deconfinement transition of the original 
gauge model (and hence of the effective spin model) is continuous, near this 
critical point, where the correlation length becomes infinite, the precise form 
of the short ranged interactions should not be important. 

d) Let us consider now the much simpler spin model with only nearest neighbour 
interactions and the same global symmetry group. If also the order- disorder 
phase transition of this model is continuous then, due to point (c) above, its 
universality class should coincide with that of the deconfinement transition. 

This last result is usually known as the "Svetitsky-Yaffe conjecture" , and has been 
confirmed in the last years by several Montecarlo simulations. It is certainly the 
most important consequence of this whole approach. 

Let us stress again that it applies only to the case in which both the deconfine- 
ment transition of the gauge model and the order- disorder transition of the nearest 
neighbour effective spin model are continuous. So it does not apply in our case 
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since in the effective actions that we shall study in the following the order-disorder 
transition turns out to be of the first order. 

There are however, in our large limit framework two other intriguing questions 
which are raised by the SY analysis. The first one refers to the difference between 
the deconfinement transition in Zat and SU(A^) LGT's. Both have the same center 
Z^r so they should be described by the same effective spin model. However we 
know very well that the two LGT's have rather different features. This differences 
becomes particularly evident in the large N limit where in one case we have the 
U(l) LGT which, for instance, in (2+1) dimensions is always confining, while in the 
other case one finds the large N deconfinement transition which will be the subject 
of the forthcoming sections. 

A second, related, problem is that, according to the SY point of view, the 
deconfined phase should correspond to the ordered phase of the effective spin model, 
but in two dimensions the Mermin-Wegner theorem [P3| tells us that the fluctuation 
are always strong enough to restore the U(l) symmetry thus forbidding the existence 
of an ordered phase and consequently of a deconfined phase in the original (2+1) 
LGT. This prediction is in clear contradiction with the analysis of the present paper 
and also with the commonly accepted scenario for non-abelian LGT's in (2+1) 
dimensions. 



2.4 Zn versus SU(iV) 

In this subsection we shall try to answer these questions and to better understand 
the difference between and SU(A^) LGT's. Let us start our discussion with 
the analysis of the Z^r models, which are simpler. In this case one easily realizes 
that the nearest neighbour effective spin model which describes the Polyakov loops 
dynamics is the Z^y symmetric "clock model". Its phase diagram in c? = 2 is well 



known ||22[. It is rather non trivial and admits, for finite N and large enough 
(3, an ordered phase in which the Zn symmetry is spontaneously broken. In fact 
the phase diagram of the Zn models (if iV > 4) is composed by three phases: 
the usual (low (3) disordered and (high /?) ordered phases and a new intermediate 
phase which is critical and has the same features of the high /3 phase of the U(l) 
model. This intermediate phase is separated from the high and low (3 phases by two 
phase transitions of the Kosterlitz-Thouless (KT) type. The Zn model in two 



dimensions is self dual, and this greatly simplifies the study of the phase diagram. 
Ordered and disordered phases are related by duality (just like in the well known 
2d Ising model) and also the two KT transitions are dual to each other [^]. The 
coupling of the dual theory is related to the inverse of the original one by a factor 
A^^. As A^ increases, the location of the first KT transition (the one which divides 
the disordered from the critical phases) remains more or less unchanged, and as 
A^ ^ oo it becomes the KT transition of the U(l) model, whose critical coupling 
Pkt is of order unity. The other KT transition and the ordered phase, which begins 
at this KT transition, due to the A^^ factor mentioned above, are pushed to infinity 
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Fig. 1: Relation between the phase diagrams of the Zjv "Clock models", the U(l) 
model and the effective model for the Polyakov loops in large N LGT at finite 
temperature. When taking the large N limit of the Z tv models, depending if the 
coupling constant is rescaled with (LGT case) or not (U(l) case), one of the 

two Kosterlitz-Thouless transitions of the Z^r models is pushed to infinity or to 
zero, respectively, and disappears for N = oo. 



as increases, and finally, in the N ^ oo limit, disappear. This agrees with the 
U(l) phase diagram discussed above, which in fact does not admit an ordered phase. 
This discussion suggests that in the large N limit the (2+1) dimensional Zjv LGT 
can never have a truly deconfined phase, and that the fluctuations of the Polyakov 
loops are always strong enough to restore the U(l) symmetry. 

In the SU(iV) case the global symmetry of the effective spin model is again Zjv, 
but the model is much less trivial. In each site instead of a single spin which can take 
N values we have a collection of N spins (the eigenvalues of the SU(A'") matrix) and 
the nearest neighbour effective action between these collections of spins is highly 
non trivial. If we label the eigenvalues with an index i then the action, which, 
according to the SY conjecture, is local in the real space, turns out to be non local 
in the index space. It is exactly this last feature which makes the difference and 
allows the existence of a high /3 phase in which the Polyakov loop expectation value 
is different from zero. Once the large A'^ limit has been taken and a master field 
configuration for the eigenvalues has been assumed, it turns out to be very difficult 
to understand how the non locality in the index space could allow to circumvent the 
Mermin-Wegner theorem. However, the above discussion on the phase diagram of 
the Zjv models allows us to gain some intuition of this phenomenon from a different 
point of view and, in particular, to see that there is no contradiction between the 
SY conjecture and the presence of a deconfined phase in the large N limit of SU{N) 
LGT's. 
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In fact, let us take a value of large, but finite and let us assume that for large 
/3 the SU(A^) LGT admits a deconfined phase. This means that all the Polyakov 
loops are "frozen" in the same direction (one of the Z^r roots of unity). If (3 is large 
enough all the eigenvalues (namely all the spins in each given site) take the same 
value and the non-locality in the index space becomes trivial. The effective action 
becomes A^^ times the action of the Zjv clock modelQ . This means that the critical 
coupling which separates the deconfined (i.e. ordered, in the language of the Z^r 
model) phase from the intermediate phase scalesQwith A^^. Thus we are exactly in 
the region where the Z^r clock model admits a broken symmetric phase and we see 
explicitly that there is no contradiction between our assumption of the presence of 
a deconfined phase and the SY dimensional reduction. Moreover we see that we can 
take smoothly in both models (LGT and spin model) the large limit and keep 
the agreement between the two phase diagrams. 

This allows us to better understand the relation between SU(A^) and Z^r LGT's 
in (2+1) dimensions. For suitable values of (5 both have the same two-dimensional 
Z^v spin model as effective action, but they are described by two very different 
regimes in the coupling space of the spin model.. Low (3 (of order unity) for the 
Ztv LGT, high (3 (of order A^^) for the SU(A^) theories. This consideration shares 
some more light on the SY conjecture and explains how is it possible that gauge 
models, which have a very different dynamics (like the Zat and SU(A^) ones) could 
be described by the same effective action. 

2.5 Character expansion in the large N hmit 

An important role in the following analysis will be played by the character expan- 
sion. Let us briefly summarize few results (for more details see Ref. 0). We shall 
particularize them to the SU(A^) case, but most of them hold for any Lie group G 
and with minor modifications also for discrete groups. The irreducible characters 
Xr{U) are the traces of the irreducible representations (labelled by r) of the group. 
They form a complete orthonormal basis for the class functions on the group. A 
function f{U) on the group is called a "class function" if it satisfies the relation: 



In particular, the characters themselves are class functions. The pure gauge action. 



^Let us stress that this is a simphfied picture and that fluctuations in the eigenvalues distribu- 
tion are present for any finite value of (3. These fluctuations are very important in the discussion 
of the phase diagram of the theory and we shall deal with them in sec. 4. But with respect to the 
present discussion they only represent a higher order correction and do not affect the validity of 
the argument. 

^ This is already apparent in the definition of the coupling (see eq. (||) ) where we have 
factorized a N"^ factor to ensure a smooth N ^ oo limit. 



f{U) = fiVUV^) 



WV e S\J{N) . 



(16) 
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The following orthogonality relations between characters hold: 

j DUxr{U)x:{U)=5r,s. (17) 

Y^drXr{UV^)=5{U,V) , (18) 
r 

where DU denotes the Haar measure (normalized to unity) on SU(A^) and dr denotes 
the dimension of the r*'^ representation. 

Besides the above orthogonality relations we shall make use of two important 
integration formulas of the characters, namely: 

j DUxr{V,U) Xs{U%) = Sr,.^^^^^ ; (19) 

I DUxriUV^U^V^) = ^X.(V^i)Xr(1^2) . (20) 
Any class function can be expanded in the basis of the characters: 

f{U)=J2Xr{U)fr , (21) 

r 

where the sum is over the set of all irreducible representations of the group, and 
the coefficients fr are given by 

fr^J DUxmfiU) . (22) 

Let us construct now the character expansion for the Wilson action. 

The Boltzmann factor associated to each plaquette in the Wilson action is (we 
neglect the index "t" or "s" of (3, which is irrelevant for the following analysis) : 

e^^^^^^^ = E^'-(/5)Xr(G'), (23) 

r 

where G, given by eq. (|]) above, denotes the ordered product of the link variables 
around the plaquette and the coefficients Fr are given by: 

/oo 
/^f/e^^^^^^^x;(t/) = E det/,^,_,+,+„(iV/3) . (24) 
n=— oo 

The rj's are a set of integers labelling the representation r and they are constrained 
by: fi > ■ ■ ■ > = 0. The indices 1 < i, j < N label the entries of the N x N 
matrix of which the determinant is taken and In{P) denotes the modified Bessel 
function of order n. 

As a consequence of the factor dr at the denominator in eq.s (^, ^) the relevant 
coefficients in the character expansion (|23|) , namely the ones that will appear in the 
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strong coupling expansions, are not the Fr themselves, but the following normalized 
coefficients: 



These coefficients have two remarkable properties, which will be important in the 
following: 

In the large N limit they have a very simple form, if the limit is taken keeping 
f3 <1 fixed: 

Fo{(3) ~ e(-) , 
FfiP) ~ ^e(-) , 

where the index / denotes the fundamental representation (whose dimension 
is A^). The above relations imply that in the large N limit 

Df{(3) = § . (26) 

Similar simplified relations hold also for higher representations. 

b) For any fixed value of N, in the large P limit the coefficients Dj. become 
equivalent to the heat kernel coefficients: 

lim Drif3) = e~W , (27) 

where Cr denotes the quadratic Casimir in the r**^ representation. As it is 
easy to see, this exponential form greatly simplifies the construction of strong 
coupling expansions. Moreover, this limit is particularly relevant for us, since 
this is exactly the situation that we have for the timelike part of the action 
when the asymmetry parameter p is sent to infinity. 



2.6 Scaling behaviour 

The ultimate test of the correctness of any lattice regularization is that, as the 
continuum limit is approached, the various dimensional quantities follow the cor- 
rect scaling behaviour. This scaling behaviour can be easily obtained by writing 
explicitly the dependence on the lattice spacing a of the relevant (dimensional) 
observables. The form of these scaling laws depends on the number of spacetime 
dimensions of the lattice. So we shall discuss separately the two cases d = 2 and 
d = 3 which are the most relevant ones for the physical applications. We shall 
concentrate in particular on the scaling behaviour of the critical temperature Tc 
(which has the dimension of a mass) since this is the simplest physical observable 
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which we can study with our techniques. We shall be interested in the scaling laws 
as functions of the coupling /5 and of the lattice size in the time direction A^^. We 
shall study for simplicity the case of a symmetric lattice Ps = Pt = P = If P 7^ 1, 
the change in the scaling laws due to the asymmetry is completely encoded in the 
equations (||) and ( p!OD and in the functions Co-(p) and Cr{p) discussed in sec. 2.1. 



2.6.1 d=2 

In the d = 2 case the coupling constant has the dimensions of a mass (see eq. 
(P)). This simplifies the analysis since in this case the coupling constant itself sets 
the overall mass scale. So near the continuum limit a physical observable, like Tc, 
with the dimensions of a mass can be written, according to the renormalization 
group equations, as a series in powers of g^. Hence in terms of the coupling (3 we 
have: 

"T.^^.^..... (28, 

The critical temperature is obtained by looking at the critical coupling P^ at which 
the deconfinement transition occurs; hence, if the lattice size in the time direction 
is Nt = we can rephrase eq. ( pS] ) as a scaling law for the behaviour of the 
critical coupling Pc as a function of A^^^. Keeping only the first term in eq. (^81) we 
find: 

PM) = a,Nt . (29) 

Thus from the simple observation that in (2+1) dimensions g'^ has dimensions of 
a mass we immediately deduce that near the continuum limit Pc must be a lin- 
ear function of Nt. If the lattice is asymmetric we can use eq.s (^ and ( p!OD to 
reconstruct the equivalent symmetric couphng and then use again eq. (|29|) . 



2.6.2 d=3 

This case is less trivial, since by looking at eq.s ^) we see that ioi d = 3 
the coupling constant g^ is dimensionless and the theory dynamically generates a 
dimensional scale Aj;, in units of which we must measure any dimensional quantity 
on the lattice. The dependence of the lattice spacing on P can be reconstructed in 
the continuum limit by using the renormalization group equations. The well known 
two loop result is: 

aAL = exp , (30) 

where bo, bi are the first two coefficients of the Callan-Symanzik equation which for 
SU(Ar) are: 

"o^i^. ^^?f^r (31) 



487r^ ' 3 Vl67r 
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Plugging eq. pO| ) into the definition of tlie critical temperature: 



Tc = ^, (32) 
aNt 



we find: 



51 



If the continuum limit is correctly taken then the ratio T^/A^ should approach for 
large enough values of P (hence, in our case, also for large values of Nt) a constant 
value. Inserting this constant into eq. (|33|) and keeping only the first perturbative 
contribution in eq. (^) we immediately see that, for d=3, j3c is a logarithmic 
function of Nt. 



3 Construction of the Effective Action 

The aim of this section is to construct an effective action for the finite temperature 
LGT in terms of Polyakov loops only. In agreement with the Svetitsky-Yaffe pro- 
gram outlined above, we shall try to integrate over the space-like variables Ui{x,t) 
in the Wilson action (^ so that the only remaining degrees of freedom will be at 
the end the Polyakov loops. The resulting effective action will live in d dimensions, 
one dimension less than the starting model. The integration over the space-like 
variables can be done in principle in two distinct steps. First one integrates over 
all the space link variables except for the ones on an arbitrarily chosen time layer. 
The result is a lattice theory with A^j = 1 in which the time-like links are the open 
Polyakov loops. This first step will be denoted for obvious reasons as "dimensional 
reduction" . The second and last step consists in the integration over the remaining 
space-like link variable, leading to an effective action with the Polyakov loops as the 
only dynamical variables. As already remarked in the introduction the integration 
over the space-like variables cannot be performed exactly and some approximation 
is needed. The approach that we shall follow consists in treating the time-like part 
St of the Wilson action as a Born term and the space-like part Sg as a perturbation. 
This means making a strong coupling expansion in jSs while treating the time-like 
part of the action exactly. The first term of this expansion, which we shall call in 
the following "zeroth order approximation" simply corresponds to neglecting the 
space-like plaquettes. In this case both steps in the integration over the space-like 
links can be performed exactly, at least as a character expansion, but the result 
is in some respect unsatisfactory. In this limit in fact the result is exactly the 



same that one would obtain with a standard Migdal-Kadanoff |^ bond-moving 
approximation and it has the same drawbacks. In fact, in 3 + 1 dimensions, it gives 
a good approximation of the whole theory only for Nt = 1 while for larger At's, 
although it still gives a good qualitative description of the phase diagram, it fails to 
predict the critical properties of the deconfinement transition, namely the scaling 
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behaviour of the critical couphng jSc as a function of Nt- This is a major problem, 
since it is only by following the correct scaling behaviour that one can finally take 
the continuum limit and extract, for instance, a reliable estimate for the deconfine- 
ment temperature. The situation is different in 2 + 1 dimensions where, at least 
to the leading order Nt-, the scaling behaviour of the "zeroth order approximation" 
and of the full theory coincide. As a consequence we may expect that the "zeroth 
order approximation" gives in this case reliable results even at large values of A^^^. 
Qualitatively the picture is the following: in 3 + 1 dimensions the statistical weight 
of the space-like plaquettes is large enough to affect the scaling properties; hence 
they can be neglected only for small values of Nt, where the critical value of j3 is 
small and the correlations between Polyakov loops induced by the space-like pla- 
quettes negligible. On the contrary in 2 -|- 1 dimensions the weight of the space-like 
plaquettes is too small to affect the scaling, and there are a good indications that 
they are not important not only for small Nt but also near the continuum limit. 

The effect of the space-like plaquettes can be taken into account perturbatively, 
order by order in (3s- In we constructed the first non trivial order in (3s for the 
SU(2) model and indeed found a better agreement with the expected scaling law in 
(i = 3. Unfortunately the effective action at this order becomes very cumbersome, 
no exact solution can be found and one has to rely on mean field estimates of the 
critical coupling. Moreover, it is easy to see that the complexity of the calculations 
further increases as higher orders are taken into account or if larger values of are 
considered. 

In the large A^ limit a completely different approach, which is more elegant 
and powerful, is available. By using suitable modifications of the Eguchi-Kawai 
techniques it is possible to perform an exact dimensional reduction, and obtain an 
action with Nt = 1 which is completely equivalent to the original Wilson action. 
Although the last step, namely the integration over the remaining space-like link 
variables, still proves to be too difficult, this is a major improvement with respect 
to previous approaches. In particular, as discussed in detail in sec. 6, the solution 
of the dimensionally reduced action in the /3s = limit gives the next to leading 
order in the scaling behaviour of (3c, which turns out to be in excellent agreement, 
in 2 + 1 dimensions, with the available Montecarlo simulations. 

This section is divided into two parts. In the first one (sec. 3.1) we review the 
construction of the "naive" effective actions, namely the ones obtained with the "ze- 
roth order approximation". In sees 3.1.1 and 3.1.2 we shall construct respectively 
the zeroth order approximation of the heat kernel action (which can be thought of 
as the p ^ oo limit of the standard Wilson action) and of the Wilson action itself. 
Then we shall obtain the corresponding scaling behaviours (sec. 3.1.3), and we shall 
discuss the reasons why they are unsatisfactory. Finally we shall describe a simpli- 
fied version of the effective action (sec. 3.1.4), which is simple enough to be solved 
exactly and at the same time accurate enough to give a good qualitative description 
of the phase diagrams of the deconfinement transition. All the material collected 
in this first part is rather old. It mainly refers to results already obtained at the 
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beginning of the eighties, even if they are discussed here in the new framework of 
our strong couphng expansion. 

On the contrary, in the second part we deal with some completely new results. 
First we shall discuss how to construct higher order terms in the jSg expansion (sec. 
3.2.1). We shall also show that, in spite of the A^^ factor in front of the space-like 
coupling, the expansion in powers of jSg is convergent order by order in the large 
limit. This is a non-trivial and important check of the consistency of the strong 
coupling expansion in (3s. Finally in sec. 3.2.2 we shall derive an exact dimensional 
reduction by using techniques typical of the Eguchi-Kawai models. 

3.1 Naive dimensional reduction {(3s = 0) and related ef- 
fective actions. 

The starting point for our considerations is the complete action defined in (^ on 
a lattice with arbitrary asymmetry parameter p. More precisely we shall consider 



e 



= n E drDM)Xr{Go.{x, t)) e^' ^^^-^^^'^^ , (34) 



x,t,i 



where a character expansion of the contributions of the time-like plaquettes has 
been performed according to eq.s (|23|) and (|25|). The naive dimensional reduction 



can be achieved by simply setting /3s = in ( p^ . In this case it is easy to integrate 
out the space link variables Ui{x,t). This can be done exactly because each link 
variable Ui{x, t) belongs only to two timelike plaquettes. In fact, consider the ladder 
of plaquettes obtained from any given timelike plaquette by moving in the time 
direction. Its contribution to the action is: 

A{x,t) =^^^Y.(^rDr{(3t)Xr{Go^{x,t))^ . (35) 

The integration over the first Nf — 1 space-like link variables on the ladder can be 
done by using eq. ([T9|), leading to the following result: 



Y^dr [Dr{(3t)f' Xr {u,{x,t = Q)P{x + i)u} {x,t = Nt)P\x)) . (36) 

■r 

The effect of this integration is to reduce the original action (which had Nt space-like 
slices in the time direction) to an effective action with only one link in the timelike 
direction and one space-like slice. The timelike plaquettes of this reduced action 
are very peculiar: the two space-like hnks coincide due to the periodic boundary 
conditions, while the two timelike links exactly correspond to the Polyakov lines 
defined in eq. (|14D , which are indeed the degrees of freedom which we are interested 
in to construct our effective action. 

The last step is now to integrate the remaining space-like links. Because of the 
periodicity in the time direction this integration is of the type given in (^). 
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By repeating the same procedure to all ladders the integration over all space-like 
link variables can be explicitly performed, leading to the following effective action 
for the Polyakov loops at /5s = 0: 

•' x,t,i I r ) 

= llJ:[Dr{/3t)f'Xr{P{S + i))Xr{P\x)). (37) 

x,i r 



Notice that up to eq. (|36|) we were still dealing with an ordinary (even if asymmetric) 
lattice gauge theory. With the last integration, the gauge theory disappears and we 
end up with a spin model in one dimension less than the original model. However 
it is still possible to recognize some remnant of the original gauge symmetry in 
the invariance of the characters under transformations of the type V —>■ W~^VW . 
Everything that has been done so far, in particular eq. (pTl), is valid for any value 
of the asymmetry parameter p. In the next subsections we shall consider the two 
extreme situations, namely p = oo and p = 1. It should be stressed that the 
effective actions obtained from (^) by setting p = oo and p = 1 are equivalent 
regularizations of the effective = theory in the continuum, but they differ away 
from the continuum limit. 

3.1.1 p — y cxD limit: heat kernel action 

The large p limit of eq. (^) can be easily obtained from the asymptotic behaviour 
of the Dj. coefficients given in (pTf), namely^ 

where we have used the fact that Nt = prit, and defined a new coupling /?hk as 

Phk{nt) = hm (39) 

This new "heat kernel" coupling is related in the continuum limit, that is for large 
Ut, to the coupling (3 of the symmetric lattice by the relation 

/5hk(^t) " Pint) + al (40) 



^(jt is a free parameter of the theory and as such is not related to Nt- It becomes a function 
of Nt if we require that as Nt changes the effective model (|37| ) describes the same physics. This 
implies for instance that the expectation value of the Polyakov loop and of other observables does 
not depend on the value of Nt- 

^Notice however that, since the space-like plaquettes have been neglected, the quantum cor- 
rections and hence the numerical values of a^'^ may be different from the ones given in sec. 2.1. 
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Away from the continuum limit, the effective action with p = oo and the one defined 
on a symmetric lattice are not analytically related, and /5hk should be regarded as 
an independent coupling. By inserting (|38D into (^) we find 

eW=o) = nEx.(P(x + i))x.(P^(^))exp (-^) , (41) 

which is the final form of the heat kernel (p oo) effective action for the Polyakov 
loop. 

3.1.2 Symmetric lattice: Wilson action 

Let us consider now the case of a symmetric lattice (p = 1), which is particularly 
important if one wants to compare the analytic results with the results of Montecarlo 
simulations. The action is obtained from ( ^7\ ) by replacing f3t with the symmetric 
couphng j3 and Nt with n^: 

^SZ,iPs=0) ^ Y[Y.[Dr{P^rXr{P{x + i))Xr{P\S)) , (42) 
x,i r 

As a consequence of eq. ( ^8|) the heat kernel and the Wilson effective actions, given 
by eq.s (^) and (^), coincide in the large /3hk and P limit provided the ratio of 
the two couplings goes to 1 in the continuum limit {nt oo). This is guaranteed 
by eq. (^) as the constant term a° is negligible in that limit. Hence, as expected, 
the heat kernel and the Wilson effective actions have the same continuum limit. 
For finite nt the two model are different and the coupling constants /3hk and /3 are 
in principle unrelated. However, in the strong coupling limit (small (3), S'p^i and 
S'pqi coincide again provided a suitable relation between the couplings of the two 
actions is established. In fact in that limit one retains in the character expansion 
of eq.s (|4l|) and (^) only the contributions of the fundamental representation and 
its conjugate, which are the leading contributions in the strong coupling limit. A 
direct comparison of the formulas, taking into account eq. (^) and the relation 
Cf = N, gives: 

P = 2e-^ + ... (/3,/5hk-^0) , (43) 

where the dots denote higher order corrections which vanish in the /5 — limit. 

3.1.3 Scaling behaviours 

This analysis and the knowledge of the explicit form of the two actions, allows us 
to obtain some general information on the scaling behaviour that we may expect 
for Tc. In the heat kernel case the result is very simple. It can be read directly from 
eq. (0) by requiring that the physics does not change if we simultaneously change 
jShk and rit. This implies a linear scaling of /5hk with nt- 

Mrit) = Jnt . (44) 
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The "renormahzed couphng" J = (3hk{nt = 1) is the relevant parameter in the 
continuum limit; it will play a major role in the following. 

In the Wilson case, due to the complicated form of the coefficients Dr{(3) the 
situation is less simple. However we can all the same deduce asymptotic scaling 



relations from the exact scaling of the heat kernel action (eq. (01)) and the rela- 



tions between (3 and (3hk (eq.s (^ , ^31) ). Let us define also in the Wilson case a 
"renormalized coupling" by setting = [3{nt = 1). We must distinguish then 
the weak coupling from the strong coupling regime. In the weak coupling the two 
actions coincide and also for the Wilson action we find a linear scaling. In fact by 
simply replacing eq. (H3) into (EO) we find: 



/3{nt) = Jnt - = (Jw + «') rit - < , (45) 

where in the last step we have anticipated the relation between J and Jw discussed 
below (eq. (0)). As already mentioned the numerical value of a° may not coincide 



with the one found by Karsch in pO[ and reported in sec. 2.1, as the contributions 
of the space-like plaquettes have not been taken into account here. The value of 
is calculated in this context in sec. 6.1 and it turns out to be close to the one 
calculated in |^ including the space-like plaquettes. In the strong coupling regime 
by using eq. (|2^) we find 

Jw = 2(^)"'. (46) 



Let us stress again that both eq. ( ^61) and the linear scaling in the weak coupling 



regime must be considered as asymptotic behaviours, while eq. ( PI) for the heat 
kernel action is exact. It is useful to write explicitly the relation between the 
normalized couplings Jw and J. They can be found immediately by setting Ut = 1 
in (|40|) in the weak coupling and from ( ^3]) in the strong coupling. In the weak 
coupling regime we have 

Jw = J-a°, (47) 

and in the strong coupling 

Jw = 2e-^ + ... . (48) 

The important lesson that we learn from this analysis is that in the continuum 
limit all these actions obey a linear scaling law. Notice that this scaling behaviour 
does not depend on the number of spatial dimensions d of the model. This is clearly 
an artifact of our approximation, and a quite disappointing one, since we know 
very well that the scaling laws do indeed depend on d. The situation is somehow 
reminiscent of what happens in the case of the mean field approximation. However, 
while in the case of mean field we know that the approximation becomes better and 
better as d increases, here the situation is exactly reverted, and the approximation 
that we make by neglecting the space-like plaquettes becomes worse and worse as d 
increases. As a matter of fact this linear scaling agrees with what we expect for the 
(2+1) dimensional gauge models, but disagrees with our expectations in the (3+1) 
dimensional case. 
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(a) 



point of 1 order transition 



line of 3 order 
transition 



1/2 



(b) 



1/(2 log 2d)~ 



1/2 



Fig. 2: (a) Phase diagram resulting from the large N limit of the effective model 
eq.^dj); (b) The same diagram expressed in terms of the heat kernel coupling J 
(related to the Wilson one by Jw = 2 exp{ — 1/(2J)}, see sec. 3.1.3, eq. (WBi)). 



3.1.4 A simplified eflFective action 



It is instructive, before analysing the rather complex model given in (0), to study 
the simplified effective action which is obtained from the Wilson one, eq. p^), by 
truncating the character expansion to the lowest order term: 



'eff 



(Jw) = E RejjwE {P{x)){PKx + i))\ 



(49) 



This model is a rather crude approximation of the original action. Due to the 
truncation, it does not even treat exactly the timelike part of the action, so it cannot 
be trusted in the weak coupling regime or in the p — > cxd limit where all the terms 
in the character expansion become important. However this model is very simple to 
study, it can be easily solved exactly and for this reason it has been very popular in 



the past years. It was extensively studied in the literature, both at finite p7| , Pq 
and in the large limit |^, with strong coupling and mean field 



approximations, and with Montecarlo simulations 



Notwithstanding its 



simplicity it turned out to be a valuable tool to understand the general features of 



the phase diagram. Its exact solution in the large N limit was derived in p9|, |30 



for any value of the space dimensions c?, leading to a phase diagram with a first 
order deconfinement transition located at Jw = It is even possible to solve 

exactly the more general model in which the Polyakov loop is also coupled to an 
external "magnetic" field /i. 



5'efr(Jw,/^) = E Re jjwE Tr(P(f))Tr(pt(f+i))| 

X \ i ) 



(50) 
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1 2 3 4 5 6 8 16 «t 

Fig. 3: Values of the critical coupling (i^ in the SU(2) theory are plotted for 
different values of the number of time-like links nt- Results obtained with Montecarlo 
simulations, which are denoted by *, are compared with those obtained in / |2^/ : A 
represents the data for fic\o, the critical coupling in the zeroth order approximation 
and □ the data for f3c\i, the critical coupling including the effect of the space-like 
plaquettes al the lowest non trivial order. 



In this framework it becomes apparent that the first order phase transition is the 



end point of a hne of third order phase transitions of the Gross- Witten type, 
located along the line 

1 — 2h 

Jw = — — , h,Jw>0. (51) 

(Jj 

As we shall see below, keeping the whole complexity of the heat kernel (or 
Wilson) effective action, the phase diagram turns out to be more rich and complex 
and also the location of the deconfinement transition changes. In particular in sec. 
5.3 we shall discuss the phase diagram of the Wilson action coupled to a "magnetic" 
field h as in eq. (0) up to the second order in the strong coupling expansion. 

3.2 Toward an exact dimensional reduction 
3.2.1 Higher orders in /3s 

We have seen in the previous subsections that the naive (zeroth order) effective 
actions fails to catch the correct scaling behaviour of the theory. Our hope is that 
we can overcome this problem by taking into account higher order contributions 



in (3s- In |^ we tested this expectation in the case of the SU(2) model in (3+1) 
dimensions and found a remarkable improvement in the scaling behaviour of the 
model by adding the terms of order This is well summarized in Fig. 3 where the 
two scaling behaviours of the zeroth order and of the (3^ order actions are compared 
with the results of the Montecarlo simulations. 
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A similar analysis in the large limit is still missing. In this subsection we shall 
give some preliminary result in this direction and in particular we shall show that, 
at least, the problem is well defined and higher order contributions do not diverge 
as N ^ oo. 



The expansion in powers of Ps of the full partition function (p4\ ) can be formally 
written as: 

g5po: ^ g5po:(/3.=0) | ^ ^ ^2^^ E E < (x, t) > 

x,t i<j 

+ In'P'sI: E E E <G,,,(fi,ti)G'.,,,(f2,t2) >+...} (52) 

Xi,ti X2,t2 ii<ji i2<j2 

where the expectation values are taken with respect to the unperturbed action, 
namely with jSg = 0. The series in brackets at the r.h.s. of eq. (|5^) can be re- 
exponentiated, leading to the following expansion for Spoi- 

Spol = Spoi^3s = 0)+N^(3sJ2J2<G^Jix,t)>, 

x,t i<j 

+ In'PsI: E E E <G,,,,(fi,tOG'.,,,(f2,t2)>c+... , (53) 

Xi,tl X2,t2 il<jl '2<j2 

where the subscript c denotes the connected part, for instance 

< G1G2 >c=< G1G2 > - <Gi >< G2 > . 

A few remarks about the general structure of this perturbative expansion are in 
order here: first it should be noticed that according to a general counting of the 
powers of in the large A^ limit we have < GiG2---Gk >c= Q( jv2fc-2 )> so that all the 
terms in the exponent at the l.h.s of (p3D are of the same order A^^. Secondly, each 
term of order involves the expectation value of k space-hke plaquettes, and a 
sum over all their possible space-time positions. If we compare different regulariza- 
tions, corresponding to different values of p, the sum over the time positions of the 
plaquettes is of order while rescales hke p"'^. So, although /3s ^ as p — > cxd, 
the effective coupling of our expansion does not vanish and it coincides with /5 + q;° 
(see eq. (0)). Indeed, as (3 is the coupling of the symmetric regularization, it is 
known from asymptotic freedom to grow logarithmically with as the continuum 
limit is approached. 

It is easy to see that the terms of order j3s in (^) vanish identically and that at 
the order j3l the only surviving contributions come from two space-like plaquettes 
in the same space position and separated by an arbitrary time interval. These type 
of contributions were explicitly calculated in the case of a gauge group SU(2) in 



Ref. |26|; for an arbitrary A^ the calculation up to order (3^ has not been done yet 



and it would involve the explicit evaluation of an integral of the type 

j DU DUU^pU\sXr{UU^)Xs{U^iU^^\) . (54) 
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This could be done by resorting to the same techniques (use of suitable Schwinger- 
Dyson equation) used in the SU(2) case |p6| . 



3.2.2 Eguchi— Kawai approach 

In the present section we shall derive an exact dimensional reduction to A''^ = 1. 
This can be achieved by using ideas similar to the ones used to reduce in the large 
N limit lattice gauge theories to twistedQ one plaquette models. This idea was first 
proposed by Eguchi and Kawai and subsequently perfected by several authors 
and it is based on the observation that in the large limit a suitably twisted lattice 
gauge theory on a lattice consisting of just one site and one link variable for each 
space-time direction generates the same set of loop equations as a theory defined 
on a large lattice, typically consisting of N^'^'^^^^'^ sites. Hence twisted one plaquette 
models can be used to describe lattice gauge theories on large lattices, by essentially 
mapping space-time degrees of freedom into internal degrees of freedom. A general 
review of the Eguchi-Kawai model, including applications to finite temperature 
lattice gauge theories, can be found in 0. Following |jl8[, rather than applying 



directly the Eguchi-Kawai method, we shall use a similar technique to reduce to 
one the size of the lattice only in the compactified time dimension. The fact that 
such exact dimensional reduction is possible is in itself an interesting result. In 
sec. B?T| we shall solve the model in the zeroth order approximation (/?<, = 0) by 



assuming as usual that the master field is invariant under translations in space. In 
this limit the final result is the same that one would obtain with a standard hot 
twisted Eguchi-Kawai model, where the dimensional reduction is done at the same 
time in all space-time direction. The present approach however has the advantage 
of being consistent in any space-time dimension, while the Eguchi-Kawai reduction 
only works for even dimensions. This will allow us to apply it in the case of (2+1) 
dimensions. 

Let us consider again the Wilson action given in eq. @. A naive prescription 
for the reduction of the degrees of freedom in the time direction would be 

V{x,t) ^ V{x) , 

U,{x,t) ^ Ui{x) . (55) 

Let us denote with 5'w(-^t = 1) the action resulting from (Q) with the substitution 
(^). It is easy to show by standard methods (see for instance ) that Sw{Nt = 1) 
leads to the same set of loop equations in the large N limit as the full theory, 
provided all loops which are closed in the reduced lattice (Nt = 1) hut correspond 
to open loops in the original lattice have vanishing expectation value. This would 
be granted by the fact that, in addition to local gauge invariance, S^^{Nt = 1) is 
endowed with the symmetry 

V{x)^e''-TrV{x) (56) 



^The twist consists in a suitable phase factor belonging to the center of SU(7V) that multiplies 
each plaquette variable in the action. 
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The trace along "open" lines is not invariant under the symmetry (^) as they do 
not contain the same number of V{x) and V^^(x) fields. So these contributions 
vanish unless the symmetry is broken. The symmetry however is actually broken 
in the weak coupling regime, that is also in the continuum limit, where V{x) is 
close to 1 (more generically to an element of Z^v) and the traces of open lines do 
not vanish. Consequently the reduction prescription must be modified, as in the 
twisted Eguchi-Kawai model, according to the formula: 

V{x,t) D{x,t)V{x)D^{x,t) , 

Ui{x,t) D{x,t)U^{x)D\x,t) , (57) 
where D{x,t) is given by 

D(f,t) = (r)S->(ro)*, (58) 

with r and Tq traceless SU(A^) matrices satisfying the 't Hooft algebra 

TTo = e'^'^ToT . (59) 

In the last equation m in an integer number to be determined. By performing in 
the Wilson action (§) the replacement (|57D and redefining the variables according 
to the substitution Ui{x) Ui{x)T and V{x) — > V"(x)ro we obtain the "reduced" 
partition function 

ll[DV{x)llDU,{x)] expiSii) , 

X i=l 

Sred = PtNY.j2^ee'"-^TT[U,{x)V{x + i)Uj{x)V\x)] + 

X i=l 

(5sN I^eTr [Ui{x)Uj{x + i)[//(f + ])U]{x)] . (60) 

X i>j 

Notice that, unlike the twisted Eguchi-Kawai model, the twists are present in (pO| ) 
only in the contributions from time-like plaquettes, as the reduction has been done 
only in the time direction. 

Consider now the loop equations for the reduced theory (|60D . We already re- 
marked that as long as the symmetry ( |5BD is unbroken the loop equations of the 
reduced theory coincide with the ones of (H). We will show now that in the twisted 
theory this is the case also in the weak coupling regime. Indeed in the extremely 
weak coupling the fields tend to their vacuum configurations, which in our twisted 
reduced theory are 

Ui{x) PN(g)lm, 

m 

Vix) ^ QN^lm, (61) 
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where Im is the m x m unit matrix and Pn and Qn_ are the usual building blocks 

m m 

for the twist eating configurations: 
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iP)ab = Sa-,i,b ; iQ)ab = 5abe'—'' , a, 6 = 1, . . . iV/m (62) 

with periodicity in the the indices a and b, namely a = ^ + 1 means a = 1. It is 
clear from the context that we have to restrict the values of so that N/m is an 
integer. With these vacuum configurations the trace of open lines is proportional 
to Tr(Qjv)*, where t is the difference between the number of V^s and y^'s in the 

m 

trace, namely the difference between the time coordinate of the initial and the final 
point of the path in original unreduced lattice. It is elementary to see from eq. ( ]62| ) 
that 

Tr(giv)* = unless t = k— (63) 

"1 m 

for integer k. In the unreduced lattice closed loops correspond to t = kNt, so the 
comparison of the two equations determines m : 

m=-. (64) 



With the above replacement eq. ( |60D is, in the large iV limit, an exact dimensional 
reduction of (0) on a (i- dimensional lattice. 



4 Phase diagram of the Effective Model 

While the exact solution of the effective model for the Polyakov loops given in eq. 
(p^) is well beyond our computational capabilities, its zeroth order approximation 
in the /3s expansion can be studied and solved analytically, within reliable approx- 
imation schemes, in both the weak coupling and the strong coupling regime. In 
particular, we consider in this section the effective action for the Polyakov loop at 
= 0, eq. (^T]), that is obtained in the p ^ oo limit (heat kernel action), with the 
aim of describing its phase diagram. 

The solution of the model (^) will rely, in both the weak and the strong coupling 
regime, on the assumption that a translational invariant master field describes the 
eigenvalue distribution of the Polyakov loop in the large N limit. In the weak 
coupling region, namely for large J = /5hk('T-t = 1), we know that the model is 
in a deconfined phase where P{x) ^ 0, and the invariant angles of the Polyakov 
loop are distributed around 6' = 0. The solution of model can then be obtained by 
retaining only quadratic fiuctuations of the eigenvalues around ^ = 0. This is done 
in subsection [4.1| , where it is shown that the model obtained from this quadratic 
approximation is a solvable model of the type known as "Kazakov-Migdal model" 
T^. The solution is a semicircular distribution of eigenvalues, centered around 



^ = 0, with a radius that increases as J decreases and acquires an imaginary part at 



a critical value of J. As first noticed by Zarembo in [|l6l , this denotes an instability 



28 



of the weak coupling solution and puts a lower bound for a phase transition. A 
deeper insight into the meaning of this phase transition, as well as a more accurate 
determination of the critical value of J, is obtained in subsection 4.1.2 following the 
remark that the contribution to the action (|4l|) from any pair of nearest neighbour 
Polyakov loops is exactly the action of QCD2 on a cylinder, of area j and with the 
holonomies at the ends of the cylinder coinciding with the open Polyakov loops. 
When the holonomies at the boundaries are given by semicircular distributions 
of eigenvalues, as in the quadratic approximation mentioned above, QCD2 on a 
cylinder can be solved exactly and it is known to have a third order phase transition 
||40[] , as Douglas and Kazakov first discovered in the case of a sphere, due to 
the condensation of instantons. In the present context this not only clarifies the 
nature of the phase transition but, as already mentioned, increases by a few percent 
the lower bound already established by Zarembo's argument. The strong coupling 
regime is studied in subsec. [4.2| . In this phase (small J) the vacuum is symmetric, 
P{x) = 0, namely it is characterised by a uniform distribution of eigenvalues. Two 
types of results can be obtained in this context. A strong-coupling expansion of 
the effective model (that was pushed up to the 4*^^ order in shows that a 
new maximum, other than the symmetric one, appears as J increases, it becomes 
competitive with the symmetric one at a critical coupling and it is energetically 
favoured for J > J^.. We have therefore, at J = Jc-, a 1^* order transition. It is 
also possible [|T^] to calculate exactly, namely at all orders in J, the mass of the 
lowest excitation in the fluctuations around the symmetric vacuum, and hence to 
calculate exactly for which value of J the symmetric vacuum becomes unstable. 
This argument determines an upper bound for the critical value of J, which, at 
least for not too high values of the number d of space dimensions, is higher than the 
lower bound determined by the weak coupling analysis, thus restricting to a rather 
narrow band the region in which the deconfinement transition must take place. 



4.1 Weak coupling expansion 

The building block of the partition function (^), that is the contribution of two 
nearest neighbour Polyakov loops, is given by a kernel of the form: 

^2{gi, 92] ^)=Y. ^r{(p) Xr{0)e-^ , (65) 

r 

where and 6 denote the invariant angles respectively of gi and g2- It is well known 
(see for instance |^) that }C2{gi, g2', A) is the partition function of QCD2 on a 
cylinder of area A and boundary holonomies gi and g2- With this notation the 
partition function (^) can be rewritten as 

Z = J llDP{x) n/C2(P(x),P(f + 0;1/J), (66) 
where J is the "renormalized" coupling J = [3t{Nt = 1) introduced in eq. (@). 
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In the large limit we can assume that the saddle point solution is translational 
invariant and, as a consequence, the action reduces to a one plaquette integral. The 
partition function is then simply given by 



DP 



/C2(P,P^;1/J) 



(67) 



where d is the number of space dimensions. Solving the model amounts to finding 
the eigenvalue distribution p{6) for the eigenvalues of the Polyakov loop P = e'^ 
that is an extreme of the free energy associated with eq. (|67D. In order to do that 
let us first rewrite eq. (|6^) explicitly as an integral over the invariant angles of P. 
This involves the explicit expression for the characters x(^) of SU(A^): 



/m _ det{e^^'^^} .n(n-i)/2 

where J^{0) is the Vandermonde determinant for a unitary matrix, 

J(e) = n2sin^, 
i<j ^ 



(68) 



(69) 



and the set of integers label the representation r of SU(iV). By using eq. (|68D 
the explicit expression of the Casimir Cr, and the Poisson summation formula 




(47rt)^/^ exp (^-ATT^nH) exp (27rin^) 



(70) 



we can rewrite (|67D in the following wayQ: 



z= mde, \j\e) 



(i-d) 



p 



(71) 



where P denotes a permutation (of signature (t{P)) of the indices. The r.h.s. of 
(fn]) depends upon a new set of integers U, which are the winding numbers of the 
eigenvalues on the unit circle^. In the weak coupling (large J) regime, we can 
assume that the invariant angles 6i are small and that the contributions of the 
winding (i.e. li ^ 0) configurations, which are exponentially depressed, can be 
neglected. We shall actually prove in the next section that in the large N limit and 

^The calculation is almost straightforward in the case of U{N) where the sum over the integers 
ri is unrestricted and the Casimir is simply given by Cr = "^iff- The details of the calculation 
in the SU(A^) case can be found in [ p6[ . 

^Similarly the integers labelling the representations can be interpreted as discretizcd mo- 
menta, and the Poisson summation formula as the corresponding discrete Fourier transform. This 



will be further clarified in the next section 4.1.2 
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above a critical value of J the vanishing of the contributions of the winding modes 
is an exact result. Hence we set /j = in eq. (^l]). Moreover, for small 9 we have 



J\9) = A\9)e-^^''^^''+''^''^ , 



(72) 



where A(^) = Y{i^j{9i — 9j) is the usual Cauchy-Vandermonde determinant. Insert- 
ing these results in eq. (|7TD we obtain 



\{d9, \a\9) 



-N[dJ-ij^\ 



n d 



(73) 



Notice that, provided J is larger than its critical value, the only approximation 
needed to go from eq. (^) to eq. (^) consists in neglecting the 0{9^) terms in 

m- 

4.1.1 Solution via the Kazakov— Migdal model 



The model given in (|73|) coincides with a Kazakov-Migdal (KM) matrix model with 
quadratic potential, which is exactly solvable. 

The KM model can be defined, in its discretized version, by the following 
partition function on a c?- dimensional hypercubic lattice: 



Zkm = J l[DU,ix) X 

X x,i 

expjj^iVTr f-V(0(f))+X:0(x)f/^(x)0(x + i)t//(x)l| , (74) 



i=l 



where (f){x) is an Hermitian N x N matrix and the link matrices Ui{x) belong to 
SU(A^). V (0) is a potential term. 

While the model (|7^) contains no self-interaction between the "gauge" fields 
Ui{x), if one integrates out the 0(x) fields the resulting effective model for the gauge 
fields (known as "induced gauge model") includes peculiar gauge self-interactions. 
The original hope that the induced gauge model could describe directly QCD in d 
dimensions was however shown to fail due to the super-confining behaviour of 
the former. 

Another way to deal with eq. ( [M]) is to carry out first the integration over 
the link matrices, to obtain an effective model for the 0(a;)'s. To do so, one utilizes 
typical matrix model techniques. First the gauge is fixed so that the 0's are diagonal: 



[X 



D(p{x) 



A\X) 



(75) 
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Then the integration over the hnks is carried out utihzing the Harish-Chandra- 

7VX;,Aa(x)ApJx+i) 



Itzykson-Zuber-Mehta integral 



-r 



A(A(f)) A{X{x + t)) 



(76) 



For a quadratic potential V (0(a;)) = |m^0^(a;), and in the large N limit, in which 
one looks for a translationally invariant master field Xa{x) Xa, the resulting model 



is: 



'KM 



ndA \A\X) 



1-d 



I P 



(77) 



Gross ]I9[ found the master field that exactly solves, for any d, this model to be 
given by a semicircular Wigner eigenvalue distribution: 

2 



P(A) 

where the radius R is given by 



- A2 



4(2rf- 1) 



(m2(c/ - 1) + c/y^m^ - 4(2c/ - 1)) ' 



(7J 



(79) 



This exact solution provides us, as we shall see in a moment, with the solution for 
the small-6' expression eq. ([73|) of our effective model. 

Let us still remark that, although unsuitable for the description of QCD in d 
dimensions, the dimensional KM model was soon argued to be related with finite- 
temperature QCD in c? + 1 dimensions P^]. Indeed the presence of "matter" fields 
0(x) in the adjoint representation of the gauge group led to conjecture that they 
could be the remnants of the components of the gauge fields in an extra compactified 
direction. This extra dimension is naturally interpreted as the time-like direction 
in a finite-temperature model. 

Going back to our effective model for the Polyakov loops, given, in the approx- 
imation of small eigenvalues, by eq. (^), and comparing it with eq. ([77D, we see 
that indeed it is a quadratic KM model^ with mass 

d-l 

(80) 



m 



2d 



6J 



The solution of the our model, which is exact as long as eq. (|73|) is exact, i.e. as 
long as quartic and higher terms in the ^'s are negligible, is therefore given by a 
Wigner distribution p{6) of radius 

2 4(2o!-l) 



J {m^{d - 1) + dJm* - A{2d - 1)) 



^The overall factor of J in the exponent in eq. (|7^) amounts just to a different normalization 
of the 6''s. 
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The radius of the distribution depends on the couphng J (a part from the overall 
factor of j) through the expression (|^) of the mass. 

Because of the dependence of rn? from the coupling J, it is easy to see that 
the argument of the square root in (^) decreases as J decreases, and it eventually 
becomes negative. Correspondingly the radius of the distribution acquires an 
imaginary part, thus making the solution of our model inconsistent. It was argued 
in [|l^] that this is a signal that the weak coupling solution becomes instable and 
that we are in presence of a phase transition. The corresponding values of J, for 
various dimensions d, are reported in the last column of Tab. I at the end of the 
following section. 

A deeper insight of this phase transition will be obtained in the next section, 
following the identification of K,2{gii 92'-, with the action of QCD2 on a cylinder. 
We just remark here that in the present calculation the instability of the weak 
coupling solution stems from the J-dependence of m?, namely from the fact that the 
Vandermonde determinant J'^[ff) for unitary matrices appears in eq. ([71|) instead of 
the usual Vandermonde determinant A^(^). So, in spite of neglecting the winding 
modes our solution "knows" that the eigenvalues live on a circle rather than on 
a line and that the excitation of the winding modes becomes eventually favourite, 
thus leading to a phase transition. 

4.1.2 Douglas— Kazakov phase transition on a cylinder 

We have already remarked that the kernel /C2(5'i, (72; defined in ( |65|) is the par- 
tition function of QCD2 on a cylinder of fixed area and boundary holonomies. The 
fact that it appears as the basic ingredient in our effective action at /^^ = is not 
really a surprise, as the space-like plaquettes are absent in ci = 1 and /3s = is exact 
in that case. So, in order to compute the partition function of QCD2 on a cylinder 
we can use a heat kernel regularization and start from eq. (|3^ ) with /3s = 0. The 
different steps of the calculation are depicted in Fig. 4, and correspond the the steps 
leading from eq. (0) with /3s = to eq. ( ^T[ ) but with the additional integration 
over all intermediate timelike links. The final integration over the remaining space- 
like link leads to the r.h.s of (|4l|) , but with the space lattice consisting just of the 
two end points of the cylinder, and with the area appearing in the exponent. This is 
exactly JC2{gi, 92', •^)- The partition function of QCD2, not just on a cylinder but on 
a generic space-time manifold is known exactly, and it has been extensively studied. 



It was soon recognised by Douglas and Kazakov that in QCD2 on a sphere a 



third order phase transition occurs at a critical value of the area. This result was 



later generalized to the case of a cylinder [^, |T5|, |10| |Tl| and the corresponding 
phase transition will therefore be called in the following "Douglas-Kazakov (DK) 
phase transition on a cylinder" . In order to study this transition let us first rewrite 
^2(91, 92', ^) by using the Poisson summation formula (|70D . The result is ^ 



>C.(,..,.A) ^ (-) e.p(^(iV=-l,)i:|^ (82) 
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QCD2 on a cylinder 



area A 



A, 



coupling P2 \ 




Polyakov loops 



coupling |32/(NjL)=P2/A 



Step 1 : 

discretization 



step3: 

^ - ^:xact integration 
over intermediate 
time-like links 



cou]3lin 



Step 2: 

exact integration 
over the space-like links 









r 









&1 



coupling P2/N|. 



Fig. 4: The partition function of the QCD2 on the cyhnder can be derived exactly 
via discretization. See the text for explanation. 



(-ir(^)5:exp 



N ^ 

[(pi - 6p{i) + 27r/i 



2A 



i=l 



where as usual P denotes a permutation of the indices. It is easy to check directly 



on ( ^21) that JC2{gi, g2', A) is the solution of the heat kernel equation on the SU(A^) 
group manifold 



- E ^^(0) - ^Ar(iV^ - 1) j /C,(0, Q,A) = ^ (8, 



uniquely determined by the condition 



3) 



B4) 



m 



6 is the invariant delta function on the group 
, redefining the Kernel by IC2{(p, 0; A) }C2 = 



where, in the notations of Ref. 
manifold. As pointed out in 
JT{4>)JT{0)K.2, eq. (^) becomes the (Euclidean) free Schroedinger equation for 
fermions on a circle, where A plays the role of (Euclidean) time. This means 
that, because of the condition (|84D , we can interpret IC2{(p,9;A) as the Euchdean 
transition amplitude for this system of fermions, from the configuration {(pi} at 
zero time to the configuration {6i} at the time A. The modular transformation 
eq. ([70D, relating eq. ( ^2] ) to eq. |6^, admits a straightforward interpretation in 
the fermionic language: the integers labelling the unitary representations in the 
character expansion ( pB| ) correspond to discretized momenta of the fermions on the 
circle, while eq. (|8^ ) gives the corresponding coordinate representation, and the 
integers in the co-root lattice are the fermion winding numbers. 
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The expressions (^,^2|) of the kernel depend upon two sets of eigenvalues, 
and {^jj.The specific eigenvalues do not tend to any limit as N goes to infinity; 
however, the corresponding eigenvalue distributions^ do have a large limit, and 
contain all the information needed to evaluate the large asymptotics of /C2. 

In particular, it is possible to write, in the large limit, the time evolution 
equation (^) as a functional differential equation for the smooth functional of the 
eigenvalue densities F[pq, pi, A], defined as }C2{(p,0;A) = exp(A^^F[po, Pi; "^]), Po 
and pi being the densities corresponding to {(/)} and {9} respectively. This goal is 
achieved by replacing partial derivatives in (|83D by derivatives with respect to the 
eigenvalue densities, 

dF d SF 



N— = — 



d(f) dx 6p{x) 



5) 



X=(f>i 



and all sums by integrals (see Ref. 
result 
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is that the time 



for the details of the calculation). The final 
area) evolution of the eigenvalue distribution 



p(x) is governed, at the leading order in by a Das-Jevicki Hamiltonian |44 



H[p{x),U{x)] 



dx p{x) 



'dU{x) 
dx 



p\x) 



(86) 



where n(a;) is the canonical momentum conjugate to p{x). In terms of the complex 
quantity 

9n(x,t) 



dx 



+ i7rp(a;,t) 



(87) 



the equation of motion becomes the Hopf equation for an ideal fluid. 



dt ~^ ^ dx^ 



with the boundary conditions 



p(x,t = 0) = vrpo(x) , p{x,t = A) = 7rpi(x) 



^9) 



The solutions to the equations ( [88|J89| ) have been studied in detail by Gross and 
Matytsin . The knowledge of these solutions allows to write the exact expression, 
in the large limit, of the free energy corresponding to the partition function ([651). 
The free energy exhibits a 3''*^ order phase transition, for a critical value Ac of the 
area [i.e. of the coupling] that we shall derive in a moment. This transition is the 
exact analogue of the 3'''^ order transition of QCD2 on the sphere, discovered by 
Douglas and Kazakov [^; actually the case of the sphere is retrieved when the 
boundary distributions both degenerate to delta functions. 



10 



If {4>i\ is a set of invariant angles, the corresponding distribution is defined as p{x) = 
Si '^(^ ~ ^i)- Notice that distributions corresponding to set of angles are periodic functions of 
x, with period 27r. 
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For our purposes it is sufficient to follow the derivation ffist given in P0[ of the 



critical value Ac- In order to find a solution to the eq.s (|8qj89| ) consider the ansatz 
corresponding to a semicircular Wigner distribution of variable radius, 



Pix,t) = — ^^r2(t) -a;2 , \x\ < r\t) 

p{x,t) = , r^{t) < |x| < vr . (90) 

By inserting this ansatz into the Hopf equation we find that it gives a solution of 
the equation provided the time dependence of the radius r of the distribution is of 
the form 

The arbitrary constants a, (3 are determined by the boundary conditions (^). For 
consistency, the initial and final distributions in (^) must be semicircular, with 
radii given by r(0) = rQ,r{A) = ri. Given tq and ri, the radius of the distribution 
is determined at any time along the cylinder by eq. (pll). However, because of the 
periodicity condition on the eigenvalue distribution, the density ( pO| , pl]) is a solution 
of the saddle-point equations (|^,^) only if r{t) < n for any t on the trajectory. 
For any given value of Tq and ri, the maximum value of r{t) increases as the area 



A increases, so the solution (^^Tj) is valid only if the area A is smaller than a 
critical value Ac, where the maximum radius equals vr and the eigenvalues fill up 
the whole circle. The critical value of A at which the transition occurs can be easily 
calculated and is given in fl^. It will be relevant in the following only the case in 
which Tq = ri = r; then 

(Acf = vr^ - vrV^ . (92) 

Notice that for r = the partition function of QCD on a sphere is retrieved. 
Consistently, the critical value (p^ ) becomes Ac = tt^, which is just the value found 
by Douglas and Kazakov [^ . 



From the previous discussion one can deduce that below the critical area Ac 
the eigenvalue distribution is confined in an interval {—a{t),a{t)) with a{t) < vr for 
any value of t. Hence the configurations in which a fraction of the eigenvalues wind 
around the circle do not contribute in the large N limit and all the integers k in eq. 
(^) can be set to zero. The topologically non trivial configuration^ are relevant 
only if at some value of t the distribution covers the whole unit circle. The role 
of instantons in inducing the Douglas-Kazakov phase transition on the sphere was 
fully investigated by Gross and Matytsin |T^, |n|]. 



A pictorial view of the DK phase transition is given in Fig. 5. 
In order to apply these results to the effective action for the Polyakov loop all 
we have to do is to remember that the role of the area is played by 1/ J and that 

These configurations are instantons in the interpretation of the eigenvalues as fermions and 
of the area as time evolution parameter. 
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Table I: Value of J below which the Wigner distribution becomes unstable, from the 
Douglas-Kazakov phase transition (third column) and from Ref. The radius 

of the distribution at the critical point is given in the second column. 



d 


rc 




j(^.c.) ]^Zarembo] 


2 


2.96 


0.321 


0.311 


3 


2.80 


0.226 


0.218 


4 


2.66 


0.192 


0.184 


oo 


0. 


I/tt^ ~ 0.101 


1/12 ~ 0.083 



the radius r(J) of the Wigner distribution of the eigenvalues of the Polyakov loop 
is given by (|HT|). Eq. ( P^ turns then into an equation for the critical value of J: 

i)' = 7r^-7rV(J)2 . (93) 

This equation can be solved for different values of the number d of space dimensions 
and the results can be compared with the ones obtained in the previous section by 
requiring the reality of r{J). The results are summarized in Tab. I. The determina- 
tion of the critical coupling obtained from the DK phase transition is consistently 
slightly higher (and hence provides a better lower limit) than the one by Zarembo. 
On the other hand the two values are very close to each other, thus indicating clearly 
that we are dealing with the same physical phenomenon, as one could argue from 
the remark at the end of the previous section, where it was noticed how the radius 
of the distribution developed an imaginary part at small J as a consequence of the 
compact and topologically non trivial support (the unit circle) of the eigenvalues of 
unitary matrices. 

The radius of the distribution at the critical DK point, given in the second col- 
umn of Tab. I, is consistently less than vr for any number of dimensions, although 
very close to it for low dimensions. The same calculation, with the critical J ob- 
tained by Zarembo, would instead give a radius larger than vr for d = 2 and d = 3. 
In spite of the small value of J we expect the weak coupling solution to be very 
reliable above the DK phase transition, particularly for large d, where the small 
value of Vc ensures that the quartic terms which have been neglected in eq. ( [75D are 
indeed small. 

To sum up, we have established that the Wigner distribution of eigenvalues be- 
comes unstable at a value of the coupling constant given by (p3|), and the transition 
is driven by the winding modes in the configuration space of the eigenvalues. This 
analogy with the Douglas-Kazakov phase transition on a sphere gives an almost 
compelling argument that the corresponding phase transition is of the third order. 
The appearance at the critical point of new classical trajectories corresponding to 
winding eigenvalues has presumably the effect of spreading the distribution, espe- 
cially at the extremes. On the other hand, for low dimensions, the transition occurs 
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Fig. 5: Mechanism of the Douglas-Kazakov transition, as it apphes to our effective 
model for the Polyakov loops 



at a radius of the eigenvalue distribution very near to vr. This probably means that 
the maximum corresponding to the broken, deconfined phase becomes unstable, 
and the distribution would collapse into the uniform one. However this is not the 
case for high dimensions, where the critical radius is small, approaching zero as d 
increases. In this case the Douglas-Kazakov phase transition is a transition from a 
classical Wigner distribution to another one, so far unknown, but still presumably 
peaked around the origin. This point will be discussed further in sec. [5.1| . 



4.2 Strong coupling expansion 

Let us turn now our attention to the strong coupling expansion of the effective model 
(p]). In this section, which is mainly based on Ref.s [|15|, [l^, we will pursue 
this expansion up to the fourth order, obtaining a satisfactory picture of the decon- 
finement transition as it is seen by approaching the critical value of the coupling Jc 
from the strong coupling side and also obtaining a rather precise determination of 
the critical value itself. 

We will moreover determine the instability of the strong coupling vacuum (i.e. 
the point where the uniform eigenvalue distributions turns from a maximum into a 
saddle point of the effective action); this instability is distinct from the first order 
deconfinement transition, as will be apparent. 



4.2.1 Notations and preliminaries 

In the strong coupling region, the vacuum of the theory is the symmetric one: 
(P) = 0; in terms of the translationally invariant large solution, it corresponds 
to a uniform eigenvalue distribution. When J grows and reaches a certain critical 
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value, a different non-symmetric vacuum, corresponding to a non-uniform eigen- 
value distribution, becomes eventually energetically favourable and the transition 
to the deconfined phase takes place. 



The partition function of the model was given in eq. (|6^) , and is repeated here 
for commodity: 



Z = J l[d9,j\e) [/C2(P,PT;1/J)]' . (94) 



We are interested in the large limit. We will therefore search for a "master field" 
solution, described by a certain eigenvalue distribution, such that it maximizes the 
free energy that appears in eq. (|9^). To proceed, we must find the strong coupling 
expansion in the large N limit of this free energy. 

As a preliminary step let us establish some notations. We are interested in the 
large N limit, so the fundamental quantity is the distribution p[9) of the eigenvalues 
of the Polyakov loop. It is convenient to expand p{6) in its Fourier modes 

1 oo 1 oo 

where an G (— 7r/2,7r/2) is the argument of p„ modulo it, and x„ coincides with 
the modulus of pn up to a sigrQ. The reality of p{9) requires p_„ = p*, while the 
normalization of p{6) to 1 in the interval (— 7r,7r) fixes po = 1. 
The inverse formula ^ 

p„ = r pie)e~'^' (96) 



shows that pn corresponds to the large limit of the loop winding n times in 
the time-like direction with the given eigenvalue distribution. In particular, p±i 
corresponds to the large A^ limit of the Polyakov loop. 

The invariance of the effective theory becomes, in the large A^ limit, a U{1) 
invariance under the shift 6^6 + 6, that is a„ — > a„ + nS. If this symmetry is 
unbroken the eigenvalue distribution is simply given by p{9) = namely = 
for n 7^ 0. In the broken phase the U{1) symmetry connects the different vacua. 
If we choose the vacuum peaked at 6' = 0, then the symmetry of the action for 
6i — > —di will force the vacuum distribution p{9) to be even in 9, thus fixing all a„ 
to zero. In this situation the eigenvalue distribution takes the form 



1 + 2 ^ x„ cos(n6') 



n=l 



(97) 



with Xn real. 



^^It is convenient for the following discussion to restrict the range for a„ and have Xn taking 
also negative values. 
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4.2.2 The integration measure 



We need to obtain a strong coupling expansion of tlie integration measure J'^iff), 
which in the large limit can be expressed as a function of the loop variables p„, 
as 



J''{e) = exp 



lim — / d9 



dip p{9) p{ip) log(l - y cos{e - if)) 



(9J 



where the double integral, which would be divergent a.t 6 = ip, has been regularized 
by the inclusion of the parameter y, and terms in the exponent suppressed by powers 
of have been neglected. Eq. (|98|) can be expressed in terms of the modes Xk by 
expanding in powers of 6 and (p and resumming the resulting expression. The result 
is 



J\9) = exp 



]imN^Co{y) + Y.Ck{y) 



X I 



k=l 



(99) 



where Co{y) is an irrelevant divergent expression and Ck{y) is given by 

k 



Ck{y) 



y 



(100) 



After removing the divergence, the limit y ^ 1 can be taken, and gives 



J'(e) = exp 



k=l 



(101) 



4.2.3 The kernel on the cylinder 

We have now to calculate the strong coupling expansion for the basic building 
block of the partition function (0), namely the QCD2 kernel on the cylinder 
/C2(P, P^; 1/./). In the following we shall consider the more general case in which 
the holonomies gi , g2 at the boundaries of the cylinder are arbitrary. 

Unlike the weak coupling regime, where we had to use the expression (^2]) for 
^2(91,92^', ^1 J) obtained by the Poisson summation formula, we use here the the 
character expansion 

/C2(^?i, 92'; 1/^) = E ^~^Xr{9i)Xr{92') ■ (102) 

R 

which provides directly a strong coupling expansion of the free energyf^ The ex- 
pansion parameter is the following function of the coupling J: 




(103) 



^•^This is the free energy of the QCD2 on the cyUnder, of course, and is just an ingredient of 
the free energy for the effective model of Polyakov loops that we are managing to build up. 
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Table II: Representations contributing to the strong coupling expansion of the 
free energy, up to the second order in . The fundamental, symmetric and an- 
tisymmetric rank 2 tensor representations must be considered together with their 
conjugate representations. 



representation 


order 


Casimir 


Character 


singlet 





1 


1 


fundamental 


1 


N 


Npi 


symm. rank 2 tensor 


2 


2N'^+2N-A 


{N'pI + Np2)/2 


N 


anti-symm. rank 2 tensor 


2 




(N^pj - Np,)/2 


N 


adjoint 


2 


2N 


N^\pi\^-1 



which vanishes exponentially as J goes to zero. The free energy F{gi, g2^', J) in the 
large limit is defined by: 

/C2(^7i,^?,-^;j) = e^^^(^-^^^^^), (104) 

where, as usual, terms suppressed by powers of in the exponent will be neglected. 
We have therefore 

F{gi,g^'; J) = 2 e-'/^'Ff{g,,g^'- J) + 4 e-'/' F^" {g,, g^'; J) + . . . . (105) 

It is found that residual dependence on J in the individual terms is polynomial in 
1/J, and hence ( |105|) is indeed a strong coupling expansion. 

The coefficients F^^ {P, P^; J) of the strong coupling expansion of the free energy 
F{P, P^; J) can be obtained in principle up to any order by using techniques which 
are summarized in Appendix A of [|15|. They involve the large A^ expansion of 



the Casimirs Cr, and the expression of the characters, in the large A^ limit, in 
terms of the Fourier coefficients pk = j^TrP^ of the eigenvalue distribution (see 
also [^5|, |12|)- In HH this program was carried out up to the 4*^^ order. The result 



is rather cumbersome and we shall not report it here, although we will use it in sec. 



4.2.4| in order to derive the critical value of the coupling with the highest possible 
precision. We just reproduce here the calculation of the free energy eq. ( |105| ) up to 
the 2'^'^ order, which illustrates all the basic points of the procedure, without being 
too involved. Let us define a representation r of SU(A^) to be of order I if the large 
A^ expansion of its Casimir is of the form Cr = IN + 0(1). It is clear from eq.s 



( |102| , |104] ) that to work out an expression of the free energy in powers of 2e 2j up 



to l^^ order, we must take into account all the representations of order up to /. The 
relevant representations up to order two are given in Tab. II, together with the 
corresponding Casimirs and the large A^ expression for their characters. 



By inserting the expressions in Tab. II into eq. (|102|) we obtain 



/C2 ~ l + 2Ar2|pi|2(l + -J^ + ...)e-^ + 
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+2-j|iVV? + iVp.P(l-^-^ + ^ + ...)e- 

+2 4|iVV?-iVp.P(l + ^-^ + ^ + ...)e- 

+ (iV>iP - l)'e~7 + 0(e-^?). (106) 

Eq. (|106|) can be exponentiated, as it is easy to check, as follows: 



/C2 ~ exp< 



r2 



In 1 f x] xi, a;?a;2 cos(2q;i — ao) \ 1 ^ 

(107) 

where we have used the notation = x^e'"'' introduced in the previous section. 
Eq. ( |107D provides the strong coupling expansion of the free energy F{P, P^; J) up 
to the second order in 2e~27. Notice that {P, P''; J), obtained from ( [L07| ) for 
A; = 1, 2, depends only on Xi and with i < k; this property is general and valid 
at all orders in the expansion. 



4.2.4 The deconfinement transition 



By using eq.s ( |101|) and (|104| , |105D we can now express the large limit of the 



partition function Z of eq. (M) as 




(108) 

It is already obvious from ( |108| ) that, for small enough J, the free energy has a 
maximum for Xk = 0, that is for a uniform distribution of eigenvalues. The first 
order approximation of ( [L08| ) can be obtained from ( [L07| ) by neglecting the second 
order terms, and reads 

-^logZ^{2e-^d-l)xl, (109) 

with all other Xk set to zero. The phase structure of (|109|) is very simple: for 
J < ^ the maximum of the free energy occurs at Xi = and we are in the 
unbroken phase. Above the critical point J = jii^^ we have instead xi = 1/2, 
which is the maximum value allowed by the positivity condition on p{6). 

A more accurate analysis of the phase diagram can be achieved by inserting 



in eq. ( |108| ) the available higher orders. As we are interested in determining the 
broken vacuum, we can set all a„ to zero, as discussed at the beginning of this 
section, and look at the maximum of the free energy regarded as a function of the 
a;„'s only. 

In looking for the maximum, it is crucial to remain within the domain where 
p{6) is positive or zero for any value of 6. This is far from trivial when we use the 
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-0.6 -0.4 -0.2 0.2 0.4 0.6 

Fig. 6: Contour plots of the free energy for: (a), J = 0.39; (b),J = 0.41; (c), 
J = 0.44. The numbers 1 and 2 mark respectively the presence of the symmet- 
ric maximum (xi = X2 = 0) and the "broken" one. The thick line encloses the 
"physical" region 



Fourier coefficients Xk as dynamical variables. The equations for such a domain can 
be obtained in principle by requiring that, for some 9q, 

p{Oo) = ^p{0)\e=e. = (110) 

and by eliminating 6q from the equations. The resulting equations for the x„'s give 
the boundaries of the physical region. 

Let us consider the second order approximation, whose corresponding free energy 
can be obtained by substituting ( |107| ) into ( |108| ). The contour plot of the free energy 



as a function of xi and X2 at various values of the coupling J is shown in Fig. 6 for 
d = 2. The physical region in the (0:1,0:2) plane can be easily determined in this 
case, and it is represented by the region inside the thick line. The straight edge 
on the right is given by the equation o;i — 0:2 = 1/2 and corresponds to density 
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distributions vanishing at = vr, whereas the the straight edge on the opposite side 
corresponds to distributions vanishing at 6^ = (Q. The plot at J = 0.39 clearly shows 
the maximum at xi = X2 = 0: the system is in the unbroken phase dominated by a 
constant distribution of eigenvalues. In the next plot, at J = 0.41, a local maximum 
has appeared at the edge of the physical region, and it is becoming competitive 
with the unbroken maximum. Notice that there is a symmetry xi —>■ —xi, so that 
a symmetric maximum appears on the other edge of the physical region. This 
symmetry is accidental, and it is removed when higher order terms are taken into 
account. In the last plot, at J = 0.44, the maximum at xi = X2 = has disappeared 
and the system is clearly in the broken phase. It is clear from this picture that 
there are three distinct regions of J: the first ranging from J = to a J^^^'"^-) where 
only the symmetric maximum exist, the second from jj^^-^'^ to some value Jjf''^'^ 
where both the symmetric and the broken vacuum are present and the third for 
J > Jlf'^'\ where only the broken vacuum survives. The first order deconfinement 
transition occurs in the middle region at the critical Jc where the value of the 
free energy in correspondence of the two maxima is the same. The point J^^^ "^ -* 
can likely be identified, when all orders are taken into account, with the critical 
coupling corresponding to the Douglas-Kazakov phase transition]^ where the weak 
coupling solution becomes unstable. The instability of the symmetric solution will 
be investigated in the next subsection, where J^^'^'^ is determined exactly, namely 
at all orders in the strong coupling expansion. 

As more orders are taken into account the explicit calculation of the free energy 
as a function of the XkS becomes more and more involved and the determination 



of the physical region more complicated. In the calculation is pushed up to 
the fourth order, where the free energy is a function of Xi,X2,X3 and x^ and the 
boundaries of the physical region are given by fourth order algebraic equations, 
together with the hj^erplanes X1 — X2+X3—X4 = 1/2 and Xi + X2 + X3 + X4^ = —1/2. 

The critical point can be determined numerically, and it is given for different di- 
mensions by the values listed in Table III, where the values of the x^s corresponding 
to the broken vacuum at the critical point are also reported. 

4.2.5 Instability of the symmetric vacuum 



In the previous section, based on |T^, we analysed the 1^* order phase transition 



separating the confined and deconfined phase. It turns out that the two phases can 



coexist as meta-stable phases, at least for not too high values of d (see sec. ^ 
Indeed the weak coupling solution is still valid for couplings J lower than Jc, down 
to the value Jj^^'^'^ determined in sec. [4.1. 2| from the Douglas-Kazakov transition. 
In this section we shall show that it is possible to determine exactly the higher limit 



^''The curved section of the boundary is part of the eUipse of equation + ~ 4x2 — and 
corresponds to eigenvalue distributions that vanish in two different points. 



^This might not be the CcLSG, cLS discussed in sec. |5.l| , for high values of d, the number of space 
dimensions. 
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Table III: Values of the critical coupling at the fourth order in the strong coupling 
expansion of the heat kernel action. In the last four columns the corresponding values 
of Xi, (i = 1 — 4) are reported. 



d 


J 


Xi 


X2 




X4 


2 


0.416 


0.41 


-0.13 


0.03 


0.07 


3 


0.282 


0.45 


-0.10 


-0.01 


0.05 


4 


0.238 


0.50 


-0.07 


-.05 


0.02 


100 


0.094 


0.64 


0.07 


-0.10 


-0.03 



jis-c.) ^Qj. ^Yie stability of the strong coupling symmetric vacuum |16, 46 , i.e. the 
value for which the configuration p = 1/2tt ceases to be a local maximum; this value 
is found to be higher than Jc. We shall give here the main points of the calculation, 



a more detailed account of it can be found in 16, 46, 17 



The goal is achieved by calculating the spectrum of excitations around the strong 
coupling symmetric, translationally invariant vacuum p{6,x) = The instability 
appears when the lowest-lying excitation becomes massless. 

As remarked several times, the configuration that dominates in the large 
limit is the solution of the classical equations of motion for the model eq. ([661). 
These equations can be derived from (|66|) by varying with respect to the eigenvalue 
distribution p{6, x) of the Polyakov loop and by taking into account the large N 
expression of the integration measure given in (|98D. After taking the derivative with 
respect to 6 , the resulting equations read: 



- P 



de'p{e',x) cot 



e-e' 



(Ill) 

where the left hand side comes from the integration measure. The dependence of 
the kernel JC2{P{x), P^(x + i); 1/J) from p{6, x) can be obtained in principle from 
the observation, discussed in sec. [4.1. 2| , that in the large A^ limit /C2(P(x), P^{x + 
i); 1/ J) is the transition amplitude, governed by the Das-Jevicki Hamiltonian (p6D, 
from an initial distribution p{9,x) to a final distribution p{9,x + i) in a "time" j. 
Hence the logarithm of the kernel /C2 is the classical action corresponding to the 
Das-Jevicki Hamiltonian, more precisely 



Ar2 



l0g/C2(P(x),Pt(f+i);l/J) 



dt / d^ 



u{e,t) 



dpie,t) 

dt 



TC 



p\9,t) 



where the field configurations are the classical trajectories, namely the solutions of 
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the Hopf equation ( ]88| ) with the boundary conditions 



(112) 



The exphcit solution of the Hopf equation with arbitrary boundary condition is in 
general not known, but we are interested in small fluctuations around the symmetric 
vacuum configuration p{6,t) = and we only need to retain in (|112|) quadratic 
terms in such fiuctuations. In this approximation the Hopf equation becomes linear, 
and hence solvable. We shall skip here the details of the calculation, which can be 
found in . Eq. ( |111| ) eventually reduces to a relation between the Fourier modes 
Pn{x) of the eigenvalue distribution p{9, x) on the sites of the lattice: 



X 



i=l 



cosh — - 
2J 



d 



— 1 n 
- — sinh — - 
d 2J 



Pn{x) + Pn{x-l) 







(113) 



which is the free field equation for the scalar excitations Pn{x) propagating on the 
lattice, provided their masses M„ are identified as follows: 



Ml = 2d (^cosh — - 1 



2(d- l)sinh 



n 
2J 



(114) 



The strong coupling vacuum becomes unstable when the first excitation becomes 
massless; setting Mi{jjf -^^^) = we obtain 



1 



X 



(s.c.) 



2iog(2d- r 



(115) 



The determination of the spectrum of excitations can be quite simply phrased also 
in the language of the previous sections [4.2.1[ - |4yT^ . Indeed it corresponds to the 
determination of the free energy ( |108| ) exactly in J, but only up to quadratic terms in 
the XnS. In particular, the determination of the strong coupling vacuum instability 
jis-c.) requires the determination of all the contributions to eq. ( |108| ) proportional 
to xl- This is within reach of the techniques developed in sec. iW^ ). and the 
result is: 



1 terms 

— - m Z — > 



2d\^e 2J — e J+e 2j — e + 
2de-^ 



d-1 



xl 



1 + e 2.J 



-1 xf 



(116) 



By requiring the above expression to vanish we get exp(— 1/2 J^^'^-'^-^) = l/{2d — 1), 



namely eq. (|115|) . 
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Fig. 7: We plot the point of instability jj^^'^'-' of the weak coupling solution, 
the point Jc^'^'^ where the strong coupling symmetric vacuum ceases to be a local 
maximum and the first-order transition point Jc in dependence of the dimensionality 
d of the lattice. It emerges that for d > 40 the weak coupling solution of eq. f p^j 
loses its validity for higher coupling than those at which the symmetric vacuum 
becomes a maximum. This points to the existence of a new phase (HI), beside 
the usual strongly coupled, confining one (I) and the weakly coupled, deconfined 
one described by the solution ( |8l| ) (II). The phases II and III are separated by a 
third-order phase transition. 



5 Extended phase diagrams 

In the previous section we studied the phase diagram of the effective model for the 
Polyakov loop, as a function of the only existing parameter, namely the coupling J. 
In this section we consider some extensions of the model involving two parameters, 
and investigate the phase diagram in their plane, which reveals some new interesting 
features. The first example we shall consider is the same model studied in the 
previous section, but regarded as a function of J and of the number d of space 
dimensions. Secondly, following Jl^ , we include the effect of external static sources 
in the adjoint representation (with coupling parameter A) and consider the phase 
diagram in the (J, A) plane. Finally we shall consider the effective action with the 
inclusion of an external "magnetic" field, as in the simplified effective action already 
studied in sec. 3.1.4, and describe how the phase diagram of Fig. 2 is modified as 
we go to higher orders in the strong coupling expansion. 
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5.1 Phase diagram in the {J,d) plane 



The results obtained in the previous section, both in the weak and strong couphng 
regime, hold for any value of the space dimensionality d. So it is possible to plot both 
j^(w.c.) marking respectively the DK phase transition and the instability 

of the symmetric vacuum, as a function of d. The plot is shown in Fig. 7. 
For low values of d (up to c? ~ 40) we have 

4™ ") <Jc< Ji"") . (117) 

This is the situation already described in the previous section: in the interval be- 
tween and jjf'^'^ both the weak coupling solution and the symmetric vacuum 
are present and a first order phase transition occurs somewhere in between. For 
d > 4:0 the situation is quite different. We have: 

Je < Ji'-'-) < Ji™-'-) . (118) 

Eq. ( |118| ) implies that if we approach the phase transition from the weak coupling 
(high J) region the semicircular solution becomes unstable before the appearing 
of the symmetric vacuum as a maximum of the effective action. This means a 
new phase, or even more than one new phase, is present for (i > 40 in the region 
/// between the two lines of Fig. 7. The features of this new, and presumably 
deconfined, phase (s) are not known, although a likely possibility is that the system 
goes from the weak coupling solution, where the density of eigenvalues vanish in 
an interval, to a distribution where it vanishes in one or in a discrete set of points. 
This phase is separated from the semicircular one by the third order transition due 
to the DK transition in QCD2. The existence of such intermediate phase(s) is also 
revealed by the circumstance that the radius of the semicircular distribution at 
the Douglas-Kazakov critical point goes to zero as d goes to infinity (see Table I), 
making the possibility of a sudden transition to a uniform distribution an extremely 
unlikely one. 



5.2 Phase diagram in presence of adjoint quarks 

Up to this point we have always been considering pure SU(A^) gauge theories. Re- 
cently in an extension of the effective model considered in section ^ has been 
considered, by taking into account the effect of a gas of non-dynamical, static 
"quarks" (that is, a gas of "electric" non-abelian external sources) transforming 
in the adjoint representation of the gauge group. The phase diagram shows a rich 



structure, similar to that of the diagram in the J, d plane discussed in sec. pA 
above. 

Static sources couple only to the "electric" link matrices V{x, t) pointing in 
the time direction. The minimal introduction of such sources is therefore that of 
a quark-antiquark pair inserted both at the same spatial location x. In presence 
of such a pair transforming in the adjoint representation, it is easy to see, by 
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using gauge invariance and the periodicity of the boundary conditions in the time 
direction, that the result of the integration in eq. (|4l|) over the space-hke hnks of 
the heat kernel action is multiplied by a factor of 

TradjP(x) = |P(x)|'-l . (119) 

The insertion of the sources is supposed to be governed by a fugacity A = e~T 
(/X being the chemical potential). Therefore, in standard fashion, one inserts k 
such pairs at locations xi,...Xk weighed with a factor of A'^/fc!, and sums over 
the positions and over k. The contribution of ( |119| ) is thus exponentiated and 
the effective model in the presence of a gas of adjoint static quarks becomes the 
following: 

Z{\) = I llDPix) n^2(P(x),P(x + i);l/J)e^S^^(l^("^)l'-^) , (120) 



which obviously reduces to the model ( |66D for A = 0. 

It is possible [|T^ to describe the phase diagram in the (J, A) plane of the above 



model. To do so, in analogy to what we did in the case A = 0, one determines 
the lines Jlf''^\X) and J^™'^^(A) at which the strong coupling and weak coupling 
solutions respectively lose their validity, and the line Jc(A) where the deconfinement 
transition takes place. It turns out that all these curves can be determined by the 
same techniques already exploited in the previous sections in the A = case. We 
will therefore presently discuss only the changes due to the new term, proportional 



to A, in the model (|I20|). 

Let us consider first the instability of the strong coupling vacuum. It is possible 
to repeat the same discussion as in sec. [4.2.5|; the extra A-term gives obviously a 



contribution, namely 

-2AIm(pi(f)e-''') , (121) 



to the eq.s of motion ( p. 1 1| ) . It is clear from ( p.21| ) that the new term affects uniquely 
the value of the mass of the 1^* excitation, that becomes: 



2(i[cosh^ - l] -2(rf- 1 + A)sinh^ (122) 



and vanishes for 

1 2rf - 1 + Arf , , 

^ = 2 log — -— . 123 

j{--)(A) ^ \-\d ^ ^ 

It is possible determine numerically the value of J at which the first order deconfine- 
ment transition occurs at fixed values of A, by adding a term d \[x\ — 1) to the free 



energy and using the method described in sec. |4.2.4 The actual values obtained 
in the case d = 2 are reported in Fig. 8. In any case, the curve of instability J'f'^'^ 
of the strong coupling vacuum represents an upper bound for the curve Jc(A) of 
deconfinement transition, as it is indicated in Fig. 8. 
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Fig. 8: The phase diagram in the A, J plane for generic dimension d is determined 
by the line J = Jc^'^'\x) where the weak coupling solution becomes unstable (1); 
the line J = j!f '^'\x) where the strong coupling symmetric vacuum turns into a 
minimum (3); the line ofl^^ order transition J = Jc(A) where the symmetric vacuum 
ceases to be energetically favourable (2). As remarked in Ref. besides the strong 
coupling confining phase (I) and the weak coupling, deconfined one (II), a new region 
appears (III), separated from (II) by a line ofy^-order transition (4). In the specific 
case d = 2 the actual numerical results for the 1^^ order transition, obtained as 
discussed in the text, are reported. The error bars estimate the uncertainty in the 
numerical procedure utilized for the calculation. 



In the weak coupling regime, the model ( |120|) can again be approximated by a 
Kazakov-Migdal model, as in the case (see sec. [4.1.1[ ) when A = 0. Indeed if we 
expand the additional term A(|TrPp — 1) up to the quadratic order in the invariant 
angles 6i of the Polyakov loop we find that eq.([75|) is replaced by: 



1 d 



E(-l)-'e 
L P 



(124) 



Eq. ( |124|) is a quadratic Kazakov-Migdal model, with A-dependent mass 

d-l\ I 



m\X) = 2ci+ 2A - 



6 



J 



(125) 



The Gross solution of this model, that describes the vacuum in the weak coupling, 
deconfined phase, is a semicircular Wigner distribution with radius r(A) determined 
by eq. (pi]), but with m replaced with m(A). As described in sec. [4.1.2| , the solution 
becomes unstable, as a consequence of the Douglas-Kazakov 3'''^ order transition in 
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QCD2, at the value that satisfies the equation 

= TT^ - 7rV(A) . (126) 




The resulting curve is plotted in the case d = 2 and sketched in the generic d case 
in Fig. 8. In analogy with the (J, d) plane, the (J, A) phase diagram reveals the 
presence of a region where neither the symmetric vacuum nor the weak coupling 
solution are stable and one or more new phases must exist. As a matter of fact, for 
A > ^, the symmetric vacuum is unstable for any value of J. 

5.3 Phase diagram in presence of a "magnetic field" 



In sec. p.l.4| we discussed a simplified exactly solvable model for the Polyakov loop. 



obtained essentially by truncating the strong coupling expansion to the first order 
and adding to the free energy a "mag netic field" term hNEx Tr(P(f) + pt(f)) I 

The phase diagram of this model consists of a line of 3^'^ order phase transition 
ending into a point of first order phase transition (see Fig. 2). The inclusion of the 
"magnetic field" to the full "zeroth order approximation" (|41| ) can be studied in 
both the weak and strong coupling regime. 

In the strong coupling expansion it corresponds to adding a term —2hxi to 
the free energy, namely to the l.h.s. of eq. ( |108D . The resulting action can be 
investigated, in the spirit of sec. [4.2.4| , by truncating for instance the expansion 



at the second order (instead of the first order as in the simplified model of p.l.4|) . 



With reference to Fig. 6 of sec. [4.2.4| , we find that the inclusion of the magnetic 
term shifts the central maximum towards the edge of the physical region even for 
J = 0: the symmetry is explicitly broken by the magnetic term. We can distinguish 
three different cases: 

a) /i<~0.15, b) ~ 0.15 </i< 1/2, c) h>l/2. 

For h < ~ 0.15 the phase transition is described by the same pattern shown in 
Fig. 6, except that the maximum "1" is not anymore located at the origin and 
it actually moves towards the edge of the physical region as J increases. In the 
simplified model corresponding to the first order truncation of the strong coupling 
expansion, this interval (/i <~ 0.15) shrinks to the point h = where the first order 
phase transition occurs. For ~ 0.15 < h < 1/2 the maximum "1" reaches the edge 
of the physical region as J increases before maximum "2" develops. The transition 
is therefore of the third order. This interval corresponds to the line of third order 
phase transition in the simplified model. Finally, for /i > 1/2 "1" is already outside 
the physical region at J = and, as in the simplified model, there is no obvious 
phase transition at this order of the strong coupling expansion|^. 

^^However, due to the non analyticity of the boundary of the physical region, a phase transition 
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In the weak coupling limit, one can reduce the model with magnetic term to 
a Kazakov-Migdal model, by expanding it in powers of 6i up to quadratic terms. 
The result is identical to eq. ( 124 ) with A replaced by h. This shows that even 
for ^ > 1/2 a phase transition of the Douglas-Kazakov type occurs, separating the 
weak coupling solution from some unknown phase. 



6 Results in the framework of the Wilson Action 
6.1 Eguchi— Kawai Model 

In sec. 3.2.2 a complete dimensional reduction was obtained by using techniques of 
the type used in the derivation of the twisted Eguchi-Kawai model. The resulting 
action (|60|) has A^^ = 1 and it is completely equivalent in the large N limit to the 
original unreduced action. By inserting the correct value of m determined in eq. 
(0) into eq. (^) we obtain the following expression for the dimensionally reduced 
effective action: 

^red = AiV^^Re e^^Tr[f/,(f)V(f + i)f/;(f)Vt(f)] + 

X i=l 

(3sNY.Y.^eTi[U,{x)Uj{x + i)Ul{x + ])U]{x)] . (127) 

X i>j 

In spite of the complete dimensional reduction, this model is still too complex 
to be solved exactly, and a solution has been obtained so far only in zeroth order 
approximation {(3s = 0) [|l^. However, even within this approximation, the solution 
will retain a non trivial dependence from Nt, resulting from the exact dimensional 
reduction. In this section we shall derive and discuss such solution. 

In order to obtain from (|127D a solvable matrix model we assume that an x 



independent master field dominates the functional integral so that the field variables 
in (|60|) can be replaced by constants. Then, setting /3s = 0, we get the following 
partition function: 

Z = JdvJIIDU, et^^NKee'^j:l^TriVU.Vfu:) ^^38) 
1=1 

For d = 3, and more generally for odd values of d (that is even space-time di- 
mensions), this result coincides with the one we would have obtained from a hot 
twisted Eguchi-Kawai model; the derivation shown here however holds for any d. 
The integrations over the unitary matrices Ui in eq. ( |128|) are all independent and 



occurs if the maximum moves from the straight to the curved section of the boundary. Numerical 
calculations indicate that this is indeed the case. Besides, the weak coupling analysis, as mentioned 
below, shows that a DK phase transition is bound to occur even for h > 1/2. 
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eigenvalue distributions \ 
become broader as P^eci^pases 




T 

2(Nrm 



Fig. 9: The weak coupling solution in the large consists of Nt distributions, 
centered at the multiples of j^. The width of the distributions grows as Pt decreases, 
until finally an instability is attained at a certain critical value Pt,c- 



the partition function can then be written as 

Z = J DY J DU gftAfRee WTr(VC/ytC7t) 

i=i 



cos(e,-0j + ^) 



where in the last hne we have gauge-rotated the matrix U to diagonahze V: 



V diag(e 



(129) 



(130) 



JT{6) is the Haar measure of SU(A^) expressed in terms of integration over the 
eigenvalues e'^' and defined in (|69|) . We want to solve the model ( |129D in the 
weak coupling limit. In the extreme weak coupling region the functional integral is 
dominated by the vacuum configuration given in eq. (|6TD , namely 



N/Nt 



(131) 



In this configuration the eigenvalues of V are organised in Nt bunches^ each com- 
posed of identical values e '^t , a = 1, . . . Nt. Quantum fluctuations around this 
vacuum configurations can be considered by parametrizing the invariant angles 6i 
according to the equation: 



6i = Or, 



2'Ka 
'N't 



(132) 



with a = 1, . . . Nt , a = 1, . . . A solution of the model given in eq. ( |128| ) can be 



obtained through the following steps: a) insert eq. (|132|) into ( p.29| ) and evaluate 
^^Notice however that the open Polyakov loop is given by V'^* 
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the integral in DU by using a saddle point method; h) in the resulting effective 
action keep only the terms quadratic in the fluctuations ^a,a] c) solve the model 
obtained in this way, which turns out to be a Kazakov-Migdal model, and find the 
the eigenvalue distribution. The details of the calculation can be found in Ref. ||18|| . 



Here we just remark about some of its most delicate points. As in the Itzykson- 
Zuber integral, the extrema of the integrand in the integral over DU are given by 
U coinciding with any element of the Weyl group, namely with any permutation of 
the eigenvalues 6'j. However in computing the integral with the saddle point method 
we only consider permutations that map each bunch of eigenvalues into the next 
one, namely: 

P ■ da,a da+l,Paia) ■ (133) 

where Pa{c() corresponds to a permutation of the eigenvalues of the bunch at 

It is easy to see that the other permutations are exponentially depressed compared 

to the ones above, typically by a factor 



exp {(3tN{l - cos^)) ~ exp L/?,iV^ ) . (134) 



for each eigenvalue mapped from a bunch a into a + 1 + s. One can assume that 
neglecting these non perturbative contributions is justified in the weak coupling 
(large jSt) regime, while below a critical value of jSt one expects a phase transition 
to occur as a result of their condensation. As we shall see, this is indeed the case, and 
we shall be able to conclude that the non perturbative contributions corresponding 
to permutations different from the ones in ( |133| ) play in the present formulation the 
same role as the instantons in the Douglas-Kazakov phase transition discussed in 
the previous sections. 

The result of the calculation, up to quadratic terms in ipa,a, is the following 
(see my- 



•' aa a ^ a,a 

( N N 

X ^ (-i)'^(^^) exp{-(AAr, - -^) Ev'a.-^a+i.P.w} • (135) 

\{P„j ^ a,a j 



The solution of this model can be obtained by assuming that the master field is 
translational invariant, namely invariant under the Z^v^ symmetry of the vacuum . 



With this assumption, which is proved a 'posteriori to be correct in [|18|, all bunches 
of eigenvalues have the same distribution and we can set: 

<^a,a ^ , (136) 



where is now of order 1 in the large Nt limit, whereas the fluctuations (y^a^^are 

Nt 



of order Notice that since the Polyakov line is P = V^^*, the angles are 
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the invariant angles of P, and the corresponding eigenvalue distribution is the 
eigenvalue distribution of P. With the position ( |136|) the partition function (|135|) 
reduces to the product of the partition functions of A^^ identical models: 



a=l 



[[ dipa [A (v^jj exp| 



Jd- 



d- 1 



12 



X 



5:(-l)'^(^)exp{nJ^<^,<^P,} 



Nt 



where we have defined n = N/Nt and done the replacement 

1 



A' 



J, 



where J is just the "normalized" coupling of the heat kernel action (see eq. 
In fact, if we consider eq. (|138|) at p = 1 we find the linear rescaling 



(137) 



138) 
I))- 



Pint) 



(139) 



in agreement with eq. (^). The comparison between eq.s ( |139| ) and (|^) also deter- 
mines the value of a^. A more complete analysis of the p dependence however can 
be obtained by writing ( |138| ) with A^^ = put and with A"^ = rit, and by eliminating 
J from the two equations. We find 



Pt{pnt)=p(^P{nt)-^^+^ . 



(140) 



This equation gives the rescaling of the coupling constant induced by varying the 
asymmetry parameter p, and it should be compared with eq. @, where Pt — ^ 
Ptipnt) and P — > Pijit). From the rescaling ( |140| ) we can read off immediately the 
function Cr{p): 

cM = -\ + ^^ + ... . (141) 
This has to be compared with the expression ( pTj) obtained by Karsch ||20 



Cr{p) = -0.2609 + — + 
Ap 



(142) 



The agreement between the two results is quite remarkable. In fact, although both 
calculations are based on evaluating one loop effects around the classical vacuum, 
this is done in the approach described here after a complete dimensional reduction 
and neglecting the effect of the space-like plaquettes, which were instead taken into 
account in Karsch's calculation. This confirms a feature, that already emerged in 
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Each eigenvalue ifa has obviously a degeneracy Nt, as it appears in each bunch. 
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the case of SU(2) namely that the corrections to Cr{p) due to the space-hke 
plaquettes are relatively small. As a consequence, and to the extent to which the 
space-like plaquettes can be safely neglected, the result ( |141| ) is independent from 
the number d of space dimensions, whereas Karsch's result was limited to ci = 3. 
The agreement between ( |141| ) and (|142| ) also represents a check of the consistency 
of the calculation, and in particular that only the fluctuations around the classical 
vacuum are relevant in the deconfined phase. 

Apart from the replacement of N with n the expression in ( |137| ) raised to the 
A^tth power coincides with the Kazakov-Migdal partition function ( [73|) obtained 
from the heat kernel action. Its solution is therefore given by the Wigner semi- 
circular distribution of eigenvalues described in section 4.3, and the instability of 
the solution occurs exactly for the same value of J as described there. The com- 
plete agreement between the results of the heat kernel action and the ones obtained 
from the Wilson action in the present section confirms that for large rit, where the 
coupling constants coincide, the two models are equivalent. This was expected, 
as discussed in sec. 3.1, and this result justifies a posteriori the assumption that 
the non perturbative contributions which have been neglected in the saddle point 
method, i.e. those involving permutations not of the form (|132|), are indeed irrele- 



vant in the weak coupling phase. It is shown in |T8| that these contributions can 
be put in one to one correspondence with the winding configurations (instantons) 
of the heat kernel model. So their condensation is expected to be responsible of the 
phase transition at the critical value of in the same way as the instantons are 
responsible for the Douglas-Kazakov phase transition. 

In terms of the original invariant angles 6i the eigenvalue distribution of our 
solution consists, as shown in Fig. 9, of a sequence of bunches, each made of 
n = N/Nt eigenvalues, at intervals of Each bunch is described by a Wigner 
semicircular distribution of radius with given by iWi). Above the phase 
transition different bunches do not communicate with each other, except for the 
canonical shift of ^ due to the twist, and the partition function is just the one 
corresponding to a single bunch raised to the power Nt (see eq. ( |137| )). The non- 
perturbative contributions, corresponding to interactions between other bunches, 
are suppressed by the exponential factor given in eq. ( |134| ). 

The picture, discussed in sec |4.1.2| , for the eigenvalues of the Polyakov loop in 
the heat kernel model is quite similar. Their distribution can be represented on an 
infinite line as a periodic distribution of period 27r, that before the phase transition 
consists of separate bunches of N eigenvalues. The contribution of a one-instanton 
configuration connecting a bunch to its neighbour is weighted by an exponential 
factor of e'XTp{(— 271"^ PtN / Nt) , as can be seen from eq. (p2|) , with A = l/J — 



ft • 

This contribution coincides with the one-instanton contribution in the one-plaquette 
model, that is obtained setting s = 1 in eq. ( p.34| ), if one takes into account that the 
number of eigenvalues in the latter case is n = N/Nt instead of A^. The contribution 
of the instantons to the classical action is then the same in the two models. This is 
a strong hint that they describe the same physics and that their condensation play 
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the same role in the phase transition. 



6.2 Strong coupling expansion 

Let us now consider the strong couphng expansion of the effective model for the 
Polyakov loops obtained from the Wilson action. We proceed by mimicking the 
treatment of the heat kernel case done in sec. ff^, and we will just sketch the main 



points, referring to for further details. 

By assuming translational invariance, we can rewrite the Wilson partition func- 
tion (0) as 

i 

whence it is evident that the integral 

X(^,Jw) = j DU exp{N J^N^e'T^Ue'^UU-'^^ (144) 

plays the same role as /C2(-P, -P^; ^/J) the heat kernel case (compare with eq. 
(|6^)). It turns out [T^ that the integral (|143|) admits a large N strong coupling 



expansion of the form: 

J(e, J) =exp(iV2f; J^Ffc(x,a)) , (145) 
^ k=i ^ 

where the functions Fk{x, a) can be shown to depend only on Xi and with i < k, 
and can be determined in principle by using Schwinger-Dyson equations. Corre- 
spondingly, the large limit of the partition function Zw becomes 

— InZw = E JwM^,<y) - -fj ■ (146) 

This is exactly the same type of expansion occurring in the heat kernel case, see 
eq. ( |108| ). The expansion parameter is now the Wilson coupling Jw instead of the 
function 2 exp[— 1/(2J)] of the heat kernel coupling. 
In particular, one finds that at the first order 

Fi{x,a) = xl = Fi{x,a;J) , (147) 

while at higher orders and F^ begin to differ]^. 

At first order therefore the phase diagram resulting from the Wilson strong 
coupling expansion is exactly the same obtained in the heat kernel case (see eq. 
( |109| ) and following text), expressed in terms of Jw instead of 2 exp[— 1/(2J)]; thus 
the critical Wilson coupling is indeed given at this order by Jw,c = 1/d. Notice that 
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The expression of Fk{x, a) up to the 4*'' order is given in |T^, eq.s (65)-(68). 
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Table IV: Values of the critical coupling Jy/^c at the fourth order in the strong 
coupling expansion of the Wilson action. In the last four columns the corresponding 
values of Xi, (i = 1 — 4) are reported. 



d 




Xi 


X2 


X3 


X4 


2 


0.601 


0.50 


0.06 


-0.03 


-0.03 


3 


0.379 


0.52 


0.03 


-0.03 


-0.02 


4 


0.275 


0.53 


0.02 


-0.03 


-0.02 


100 


0.010 


0.49 


0.004 


0.004 


0.005 



this is, as it was to be expected, the result that is found within the simple model 
reviewed in sec. 3.5, a model that was indeed obtained from the truncation to the 
first order of the Wilson action. 

Of course the deconfinement transitions observed in the Wilson and the heat 
kernel regularization schemes correspond to the same physical phenomenon. From 
the above discussion it follows that this is the case if 

Jw = 2exp(^-i^) + ... , (148) 

the extra terms in the r.h.s. being the effect of higher order corrections. This is 
just the relation between the two couplings given in P8| ) and valid in the strong 
coupling regime. This relation was derived in sec. 3.3 by comparing the strong 
coupling behaviour of the coefficients in the character expansions of the heat kernel 
and Wilson actions. The analysis of the deconfinement transition, by using the 
strong coupling expansion of the Wilson action, can be carried on in exactly the 
same way as in sec. 4.2.4, and the results up to the fourth order are summarized in 
Table IV. The results of Table IV all refer to nt = 1. It is possible in principle to 
apply the same method to extract from the dimensionally reduced twisted action 
of eq. ( |128| ) with p = 1 a strong coupling expansion for any value of Uf. However, 
as the Polyakov loop coincides in ( |12^ ) with TrV^"% namely with x„„ higher and 
higher orders in the strong coupling expansion are needed to proceed to higher 
values of rif. It is more convenient to use the values of Table IV as an input and 
use the scaling law ( p.39| ) to extract the dependence from rit. 



7 Comparison with results from Montecarlo sim- 
ulations 

When comparing our results with those of Montecarlo simulations we must expect 
two kinds of systematic deviations. 
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a) The first one is due to the large N approximation. This is however a rather 
small deviation. In fact, the lesson that we learn by looking at the Montecarlo 
simulations is that the large N limit results are reached for rather small values 
of N. For instance recent calculations on the glue-ball spectrum in (2+1) 
dimensions for low values of have shown that already for N > 3 some 



mass ratios are well described with the leading term in the large N limit ||47|| . 
Thus we can consider our large N results as rather good approximations of 
the = 3 case in which we are interested. At the same time rather large 
corrections are certainly present in the N = 2 case. Obviously, it would be 
highly desirable to have results from simulations for larger values of A^, so 
as to study in detail the approach to the large A^ limit and test the analytic 
calculations reviewed here in a more stringent way. But this is not an easy 
task. Large A^ simulations have been performed only in the framework of the 
Eguchi-Kawai reduction scheme and are plagued by the presence of a bulk 
phase transition which shadows the deconfinement one. 

b) The second, more serious, problem is that, in order to solve the model, we had 
to neglect the space-like plaquettes in the original action. We have already 
discussed the consequences that this approximation has on the scaling laws 
in sees 3.1.3 and |6]^: only for d = 2 we find the correct scaling behaviour. 
The experience gained with the SU(2) p6| model however suggests that in 
the rit = I case the corrections due to the space-like plaquettes are very 
small and that their effect is essentially negligible. Even if this case is not 
interesting from a phenomenological point of view, since rit = 1 is too small 
to give informations on the continuum limit, it becomes very interesting as a 



test of our method. For this reason we shall devote sec. to a discussion 
of this comparison. For larger values of these corrections are not any 
more negligible and their importance increases as rit increases. Notice that, 
as mentioned above, the corrections due to the inclusion of the space-like 
plaquettes have a different relevance in (2+1) and (3+1) dimensions. In (2+1) 
dimensions they do not affect the scaling behaviour, and in fact we shall be 
able to successfully compare not just the leading, but the next to leading term 



in the scaling law ( |139[ ) with the results of Montecarlo simulations. On the 



contrary in (3+1) dimensions they completely change the scaling law. For 
this reason we shall discuss separately the (2 + 1) and the (3 + 1) dimensional 
cases in the following subsections. 

7.1 rii = 1 in (3+1) dimensions 

In the Ut = 1 case we have only two results from Montecarlo simulations; they 
are however very interesting and carry a lot of information. The first one is for 
the SU(2) model It is very precise since it was obtained by using an original 



and very powerful non local algorithm (see |^ for the details). The value of 
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the critical temperature in our units is Jw,c = I3c/N'^ = 0.2185(1). The second 
results is for = 30, hence in a situation in which the large limit is a very 
good approximation. The critical coupling is Jw,c ~ 0.38. It was obtained in [ES 



within the framework of the Eguchi-Kawai reduction scheme. These results must 
be compared with what we have found in the case of the Wilson action. In this 
case our best estimate for the critical coupling in (3 + 1) dimensions, obtained by 
means of the strong coupling expansion, is Jw,c = 0.378 as shown in Table IV. The 
agreement with the = 30 value is indeed impressive. This is an interesting result 
and, compared with the naive results Jw,c = l/d = 0.333 which one would obtain 
by using the simplified action of sec. 3.5, it shows the relevance of the higher order 
terms in the character expansion of our action. The rather large gap between the 
value of Jw,c for N = 2 and = 30 can be understood already at the level of the 
lowest order in the strong coupling expansion where, for N = 2 only, the reality of 
the fundamental representation leads to a critical coupling Jw„c = It would be 
interesting to have the Montecarlo simulations for some intermediate value of A^, to 
see how the large A^ limit is approached. Analytically this result could be obtained 
within the strong coupling expansion scheme developed for the large A^ limit. A 
different technique has been applied both for N = 2 [EH] and for higher values of A^ 



5% consisting in an improved mean field Bethe approximation that determines the 
asymptotic expansion of Jw„c in powers of l/d. The results obtained so far can be 
summarized as follows: 

iV>3 =l+-L + 0(i). (149) 

Notice that for A^ larger than 3 the expansion is the same up to terms of order ^, 
which means that already at A^ = 4 we are quite close to the large A^ limit. If one 
substitutes (i = 3 in ( |149| ) keeping all known terms one obtains the values 0.2045, 
0.306 and 0.389 for A^ = 2, A^ = 3 and A^ > 3 respectively. They significantly 
improve the lowest order results and they agree, within the approximation one 
would expect from the truncation, with the Montecarlo results. 

7.2 2+1 dimensions 

Very precise Montecarlo estimates of the deconfinement temperature in (2+1) di- 
mensions exist for the N = 2 and A^ = 3 models in the range 2 < < 6. These 



can be found in [g^, |51|, ^ and are reported in Tab. V. All these simulations were 
made with the standard Wilson action and at p = 1. 

Looking at the data one can see that the expected linear dependence on nt is 
very well fulfilled. This allows us to use the variable = ^, which is the natural 
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one to take the continuum limit (see Tab. VI). Once the hnear term is factored 
out we can see in Tab. VI a residual dependence on nt which is due to the higher 
order terms of eq. (pSf). These corrections are so small that one can safely use a 
simple one parameter fit, keeping only the leading 0(1/ (3) correction, to extract 
reliable estimates for the continuum limit values of the critical temperature. These 
are reported in the last row of Tab. VI 

As discussed in the previous sections, in the framework of the EK reduction 
scheme we can predict both the leading and sub-leading scaling behaviour. Re- 
markably enough, not only the functional form of these two terms, but also the 
numerical estimates of the correction that we find are in very good agreement with 
the simulations. They are reported in the last column of Tab. VI. They have been 
obtained by taking for /?c(^t = 1) = Jw,c our best strong coupling estimate given 
in Tab. IV and by using the scaling law ( |139| ) to extrapolate the result to arbitrary 
nt- The result is: 

- - 0.351 + . (150) 



Ng"^ nt Ant 

This analytical prediction has to be compared with the following best fits, for the 
N = 2 and = 3 cases, obtained from the Montecarlo data of Tab. VI by 
assuming a linear dependence of the critical coupling from 1/nt. Both fits have a 
high confidence level and give: 

JV = 2 Mm) ^ 0.380(3) + "°^'"> . 

* nt 

AT = 3 ^^ = 0.366(2) + HUM . (151) 

nt 

The agreement of the analytic results (|150|) with the best fit (|151|) of the Mon- 
tecarlo simulations is quite remarkable. In fact the trend shown in the N = 2 and 
A^ = 3 simulations is perfectly consistent with the theoretical large A^ limit not 
only in the leading term but also in the sub-leading one. Although it is premature 
to draw any conclusion, one can say that the results are at least consistent so far 
with the space-like plaquettes being negligible in 2 + 1 dimensions as far as the 
determination of the deconfinement transition is concerned. Notice that in c? = 2 
there is no simulation for = 1 to make a direct comparison as we did in the (3+1) 



case above. In Q however we obtained a rather reliable estimate of (3c iii the 
SU(2) case. This value is reported in the first line of Tab. VI. 

7.3 nf > 1 in (3+1) dimensions 



Very precise data on SU(2) and SU(3) in the range nj = 2 — 16 can be found in |]54 



Some results for higher values of A^, obtained within the context of the EK model 



can be found in [^, We have collected in Tab. VII these results for the lowest 
values of nt- 
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Table V: The critical coupling (3c as a function of the lattice size in the t direction, 
nt, in the (2+1) dimensional SU(2) and SU(3) LGT, taken from /Q/ and /Q/. 





= 2 


= 3 


2 


~ 0.866 


0.906(2) 


3 


~ 1.251 




4 


1.630(8) 


1.638(6) 


6 


2.388(10) 


2.371(17) 



Table VI: The critical temperature = in (2 + 1) dimensions, as a function 
of nt and N. In the first line of the N = 2 column is reported an analytic prediction 
obtained in for nt = 1, and in the last line the extrapolations to the continuum 
limit. The prediction of the analysis reviewed in this work is displayed in the last 
column. 



nt 


N = 2 


N = 3 


= oo 


1 


~ 0.459 




0.601 


2 


~ 0.433 


0.4531(8) 


0.476 


3 


~ 0.417 




0.434 


4 


0.4075(19) 


0.4094(14) 


0.413 


6 


0.3979(16) 


0.3952(27) 


0.393 


oo 


0.380(3) 


0.366(2) 


0.351 



Table VII: The critical coupling f^c/N"^ , in (3 + 1) dimensions, as a function of ut 
and N. In the last column our predictions, according to the scaling law eq.(p^ and 
PJN^ = 0.378 fornt = 1. 



nt 


iV = 2 


A^ = 3 iV = 24 7V = 54 


TV = 81 


TV = 96 N = oo 


2 


0.4700(8) 


~ 0.568 - 0.70 




0.506 


3 


0.5442(8) 


0.6167(11) - 0.70 




0.634 


4 


0.5746(2) 


0.63250(2) 




~ 0.70 0.762 


5 


0.5932(11) 






0.890 


6 


0.6066(8) 


0.65490(6) 


0.705(5) 


1.018 



As mentioned above in this case we do not obtain in our approach the correct 
scahng law for the critical temperature as a function of rit for which the inclusion of 
the space-like plaquette contribution is mandatory. Hence we can trust our results 
only for low values of We have already seen before that in the = 1 case our 
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prediction is indeed successful, we shall now extend the comparison to higher values 
of rit. To obtain our best estimate for the critical temperature in this case we use 



again the scaling relation eq. (|139|) and take our rit = 1 value ( Jc = 0.378) as input. 



These estimates are reported in the last column of Tab. VII. 
A few comments are in order at this point: 

a) It is interesting to notice that the critical couplings for small values of nt 
have a dependence on rtt which is not too different from ours. In fact, even if 



the scaling law eq. (|139| ) is definitely different from the scaling behaviour eq. 
(|33D, for low values of rit the Montecarlo estimates of the critical couplings 
are still far from the asymptotic scaling region and do not follow the scaling 
behaviour of eq. (|3^) . 



h) At fixed rij, as increases also the critical temperature systematically in- 
creases, and our large estimates seem to be the upper limit of this be- 
haviour. 

c) For large values of A^ all the Montecarlo result cluster around the value ~ 0.70. 
This is an artifact of the EK approximation. This value is due to the presence 
of a bulk transition which shadows the real deconfinement transition. 

As we have seen, our results, obtained analytically in the large A^ limit, show a 
more than qualitative agreement with Montecarlo simulations performed on the 
full theory. We think that a careful analysis of higher order contributions in /3s, 
may allow in future further improvements in the analytic determination of the 
deconfinement critical temperature. 
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